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I.  EXECUTIVE  SUMMARY 


Motivated  by  reports  of  accidents  in  high- altitude  aircraft  flights  due  to  unforecast  strato¬ 
spheric  turbulence  caused  by  gravity  waves  over  mountains,  this  research  program  aimed 
to  improve  the  theoretical  understanding  of  various  physical  mechanisms  relevant  to  the 
generation  of  stratospheric  gravity  waves  by  orography.  Specihcally,  the  main  objectives 
were  to  study:  the  effects  of  the  Earth’s  rotation  on  nonlinear  mountain  inertia-gravity 
waves;  the  generation  of  gravity  waves  by  transient  wind  over  mountains;  the  interaction  of 
mountain  gravity  waves  with  the  tropopause  and  possible  amplihcation  of  these  waves  in  the 
stratosphere;  and  possible  transfer  of  energy  from  shear-flow  instabilities  to  low-frequency 
inertia-gravity  waves  via  resonant  nonlinear  interactions. 

Towards  reaching  these  goals,  the  following  problems  have  been  studied: 

(i)  Effect  of  the  Earth’s  rotation  on  nonlinear  mountain  waves.  Eormally,  the  effect  of 
rotation  on  the  generation  of  mountain  waves  by  wind  is  controlled  by  the  Rossby  number, 
Ro  =  U /  {Lf),  where  U  denotes  the  wind  speed,  L  is  the  mountain  length  and  /  the  Coriolis 
parameter.  Rotation  is  expected  to  be  relatively  unimportant  when  Ro  >  0(1),  and  this 
regime  is  realized  if  the  mountain  is  not  too  long,  L  <  50  km  typically.  However,  even  when 
Ro  3>  1,  rotation  does  affect  the  low-frequency  part  of  the  gravity  wave  spectrum,  suggest¬ 
ing  that  weak  rotation  effects  could  become  important  at  long  distances  from  the  mountain. 
This,  in  turn,  raises  the  question:  under  what  conditions  can  rotation  be  safely  neglected  in 
analyzing/modeling  mountain  waves? 

We  have  addressed  this  issue  by  constructing  an  asymptotic  model.  It  turns  out  that  ro¬ 
tation  behaves  as  a  ‘singular  perturbation’:  the  effects  of  rotation  become  equally  important 
to  the  effects  of  stratification  far  from  the  mountain,  regardless  of  how  large  Ro  is.  As  a 
result,  the  response  far  from  the  mountain  is  dominated  by  inertia-gravity  waves  that  would 
be  absent  had  rotation  been  ignored.  Hence,  the  wave  signature  of  the  mountain  extends 
much  further  than  would  be  expected.  This  piece  of  work  forms  part  of  the  doctoral  thesis 
of  B.  Druecke  (in  preparation).  Eor  technical  details,  see  §1  in  Sec.  H  below. 

(ii)  Combined  effect  of  rotation  and  the  tropopause.  The  asymptotic  model  above  was 
generalized  to  account  for  the  combined  effect  of  the  Earth’s  rotation  and  the  tropopause. 
In  the  presence  of  the  tropopause,  the  induced  wave  disturbance  over  the  mountain  can  be 
dramatically  increased  when  the  tropopause  is  at  certain  heights  (depending  on  the  wind 
speed  and  the  stratihcation).  Under  such  ‘tuned’  conditions,  nonlinear  wave  interactions  are 
particularly  strong  and,  combined  with  rotation,  drive  inertia-gravity  waves  that  propagate 
far  downstream  of  the  mountain.  Our  theoretical  results  indicate  that  this  generation  mech¬ 
anism  is  robust  for  mesoscale  mountains,  and  rotation  effects  can  be  far  more  important 
than  previously  thought,  especially  when  tropopause  tuning  is  taken  into  account.  A  nu¬ 
merical  code  was  also  written  for  solving  the  Euler  equations  in  the  presence  of  stratihcation 
and  rotation.  Fully  numerical  simulations  of  unsteady  responses  for  various  heights  of  the 
tropopause  lend  support  to  the  theoretical  predictions.  This  piece  of  work  is  also  part  of  B. 
Druecke’s  doctoral  thesis  (in  preparation).  Technical  details  can  be  found  in  §2  of  Sec.  H. 
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(iii)  Effect  of  unsteady  wind  on  nonlinear  mountain  waves.  Previous  work  on  mountain 
gravity  waves  invariably  assumes  steady  wind;  apart  from  a  few  studies  of  the  effect  of  a  time- 
harmonic  wind  component  on  the  linear  response,  the  possible  signihcance  of  wind  variations 
remains  unexplored  theoretically.  Using  an  earlier  study  by  Kantzios  and  Akylas  (Proc.  Roy. 
Soc.  bond.  A  440,  1993)  as  the  starting  point,  an  asymptotic  model  was  developed  for  the 
evolution  of  the  nonlinear  response  in  the  presence  of  slowly  varying  wind  (timescale  of  the 
order  of  hours).  In  the  case  of  time-harmonic  wind  variation,  the  model  reveals  that  the 
appropriate  perturbation  parameter  is  the  relative  amplitude  of  the  oscillatory  wind  compo¬ 
nent  divided  by  the  wind  frequency  (made  dimensionless  by  the  buoyancy  frequency).  As 
a  result,  a  small-amplitude  low-frequency  oscillatory  wind  component  can  play  a  signihcant 
role,  and  the  response  may  be  quite  different  from  that  predicted  by  steady-state  theory. 
We  have  also  developed  a  fully  numerical  code  to  study  the  effects  of  unsteady  wind  and 
validate  the  theoretical  predictions.  This  work  is  part  of  B.  Druecke’s  doctoral  thesis  (under 
preparation).  Technical  details  are  given  in  §3  of  Sec.  II  below. 

(iv)  Energy  transfer  to  low-frequency  gravity  waves  via  resonant  nonlinear  interactions. 
We  have  examined  the  radiation  of  low-frequency  gravity  waves  by  weakly  nonlinear  wavepack- 
ets  propagating  in  a  stratihed  fluid  of  variable  buoyancy  frequency.  When  wave  trapping 
is  possible  as,  for  example,  in  the  case  of  a  density  inversion,  we  hnd  that  the  mean  flow 
induced  by  nonlinear  interactions  can  give  rise  to  radiating  gravity  waves.  This  mechanism 
is  analogous  to  the  transfer  of  energy  from  instability  wavepackets,  generated  by  shear-flow 
instability,  to  low-frequency  radiating  gravity  waves.  This  is  work  still  in  progress. 

The  main  results  of  the  present  study  are  highlighted  and  discussed  in  connection  with 
atmospheric  gravity  wave  generation,  in  §4  of  Sec.  II. 

•  Personnel  supported:  Professor  T.  R.  Akylas,  B.  Druecke  (graduate  research  assistant) 
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II.  TECHNICAL  DESCRIPTION 


1  Rotation  Effects  on  Nonlinear  Stratified  Flow  Over  Topography 

In  the  atmosphere,  the  typical  value  of  the  Brunt-Vaisala  (buoyancy)  frequency  N  ~ 
10“^  rad/s  and  the  wind  speed  17  ~  10  —  30  m/s,  so  the  characteristic  vertical  length- 
scale  of  mountain  waves,  17/A^  ~  1  —  3  km,  is  much  shorter  than  the  typical  length  of  the 
topography,  L  ~  10  —  100  km;  this  justihes  using  the  hydrostatic  approximation,  taking 
dispersive  effects  to  be  weak.  The  signihcance  of  rotation  effects  is  controlled  by  the  Rossby 
number  Ro  =  U/  {Lf),  f  being  the  Coriolis  parameter,  and  f/N  ~  10“^  in  the  atmosphere. 
On  scaling  grounds,  then,  rotation  would  be  expected  to  be  relatively  unimportant  when 
Ro  ^  1  and  this  regime  is  realized  if  the  topography  is  not  too  long,  L  <  50  km.  Here  we 
make  a  systematic  study  of  gravity  waves  over  topography  in  the  limit  Ro  3>  1,  which  reveals 
that  weak  rotation  can  have  a  signihcant  effect  in  the  far-held  response.  For  simplicity,  we 
hrst  consider  the  case  of  N  constant.  The  effects  of  the  tropopause,  which  turn  out  to  be 
crucial,  are  included  in  §2. 

1.1  Governing  Equations 

The  governing  equations  for  incompressible,  inviscid,  stratihed  how  in  a  rotating  coordinate 
system  are 


V-u  =  0 


Rp 

Dt 


=  0 


where  /  =  212  sin  0  is  the  Coriolis  parameter,  12  =  7.2722  rad/sec  being  the  angular  speed 
of  the  earth’s  rotation  and  6  the  latitude.  The  Coriolis  parameter  ranges  in  value  from  /  = 
0  sec“^  at  the  equator  where  rotation  does  not  locally  ahect  hows  to  /  =  ±1.4544  x  10“^  sec“^ 
at  the  poles  where  local  rotation  ehects  are  maximum.  The  assumption  that  the  thickness 
of  the  huid  is  relatively  small  compared  with  the  horizontal  length  scales  has  been  made  fol¬ 
lowing  Kundu  and  Cohen  [12].  The  coordinate  system  is  chosen  such  that  y  is  the  outward¬ 
pointing  normal  to  the  earth’s  surface  (vertical),  x  points  in  the  eastward  direction  and 
southward,  as  depicted  in  Fig.  1. 


Assuming  the  how  is  two-dimensional  such  that  there  is  no  variation  in  the  direction, 
the  velocity  components  in  the  streamwise  and  vertical  directions  can  be  written  using  the 
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Figure  1:  Local  coordinate  system  xyz  rotating  about  polar  axis  with  frequency 


streamfunction,  as 

\i{x,y,t)  =  {u,v,w)  = 

The  governing  equations  then  become 

pt  +  J  (p,T)  =  0 
p  {^yt  +  J  (Ty,  T)  +  fw}  =  -p^ 
p{-'^xt  -  J  =  -Py-  9P 

p{Wt+J  {w,  T)  -  f^y}  =  -P^ 


(2) 


(3) 


where  J  (0,  p)  =  4>xVy  ~  4>yVx  is  the  Jacobian.  For  inviscid  flow,  only  the  kinematic  boundary 
condition  can  be  imposed  at  the  surface  of  the  topography.  For  two-dimensional  topography 
given  hy  y  =  h  {x)  as  depicted  in  Fig.  2,  the  kinematic  boundary  condition  is 

(JT 

u  ■  n  =  0  on  y  =  h{x)  ^  —  =  0  on  y  =  h{x)  T  =  const  on  y  =  h  {x)  (4) 

(JjJb 

The  assumption  that  disturbances  do  not  propagate  far  upstream  is  made.  Therefore,  far  up¬ 
stream  as  a;  — — cx),  the  wind  velocity  is  uniform  giving  T  — Uoy  and  there  is  no  transverse 

dp 

velocity,  w  ^  0.  The  density  held  is  unperturbed,  p  ^  p{y)  and  hydrostatic,  —  =  —gp. 

dy 

Finally,  the  pressure  gradient  in  the  transverse  direction  is  geostrophic,  giving  p^  =  pfU^. 
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Figure  2:  Physical  system  of  stratified  flow  over  topography  in  the  presence  of  rotation. 


1.1.1  Nondimensionalization 

We  nondimensionalize  the  variables  as  follows 

X  =  Lx',  y  =  ^y',  t  = 

iVo  Uo 

U 

'h  =  w  =  fLw',  p  =  pop',  p  =  gpo^p' 

l\o  l\o 

where  L  is  the  characteristic  width  of  the  mountain,  Uq  the  characteristic  wind  velocity,  Nq 
the  characteristic  Brunt-Vaisala  frequency  and  po  the  characteristic  density.  Substituting 
these  into  the  governing  equations  of  (3)  and  dropping  the  primes  gives 

Pt  +  J  (p,d')  =  0 


Pp  +  J  (d'y,  d')  +  =  -Px 

PP^P  +  J  ('ha;,  \h)}  =  Py  +  P 

p{Wt  +  J  {w,  'h)  -  d'y}  =  -p 

where  the  dimensionless  parameters  are  the  Boussinesq  parameter,  P  = 


(6) 


UqNq 
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,  the  longwave 


parameter  p  = 


and  the  Rossby  number  Ro  =  The  Boussinesq  parameter  is  the 
NqL  j  L 


ratio  of  the  vertical  length  scale  of  fluid  motion  to  the  vertical  length  scale  over  which  den¬ 
sity  changes,  and  the  limit  as  0  gives  the  Boussinesq  approximation.  The  longwave 
parameter  is  the  ratio  of  vertical  to  horizontal  length  scales,  and  the  limit  as  p  — 0  gives 
the  hydrostatic  approximation.  The  Rossby  number  is  the  ratio  of  the  Coriolis  time  scale 
to  the  advective  time  scale.  Thus  in  the  limit  as  Ro  oo,  rotational  effects  are  relatively 
unimportant.  The  limit  where  Ro  =  1  is  the  limit  where  both  rotation  and  advection  are  of 
equal  importance,  and  as  Ro  0,  rotation  dominates. 


Far  upstream  the  flow  is  undisturbed,  and  the  nondimensionalized  flow  is 


y,  w  ^0,  p  =  -ppN  {x 


-oo 


(7) 
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The  nondimensionalized  boundary  condition  on  the  topography  becomes 


T  =  const  on  y  =  eh  {x) 


(8) 


where  e  =  — — ^  is  a  parameter  determining  the  importance  of  nonlinearity  and  ho  is  the 
Uq 

characteristic  height  of  the  topography.  The  nondimensional  topographic  prohle  used  in  the 
analyses  summarized  below  is  either  the  Gaussian  prohle  of  functional  form 


h{x)  =  e 

or  the  algebraic  (Witch  of  Agnesi)  prohle 


(9) 


h  (x)  = 


1  + 


(10) 


The  how  is  assumed  to  be  hydrostatic,  implying  /x  — 0,  and  nearly  steady.  We  dehne  a  slow 
time 

T  =  z/  1  (11) 

Then  the  time  derivatives  in  the  governing  equations  can  be  rescaled  and  the  hydrostatic 
approximation  imposed  to  give 

J(p,  T)  =  -z/Vr 


(3p 


=  -Px 


Py  =  -P 


(12) 


p  |z/^wt  +  j  (w,  th)  -  =  -p 

Diherentiating  the  second  equation  with  respect  to  y  and  eliminating  pressure  in  favor  of 
density  using  the  third  equation  gives 

J(p,  T)  =  -z/Vr 


(^{p 


z/^T^y  +  J  (T^,  T)  + 


Rd^ 


w 


Px 


X  J  y 


Py  =  -P 

J  (tc,  T)  =  -  (l  +  v‘^wt) 


(13) 
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1.1.2  Integration  of  Density  Eqnation 

Making  the  assumption  that  the  streamlines  never  become  vertical,  '^y  7^  0,  throughout  the 
domain,  one  can  make  the  coordinate  transformation 

{x,y,T)  ^  .  (14) 

Then  substituting  this  into  the  density  equation  above  gives 

Px\q,  =  (15) 

Dividing  through  by  and  integrating  with  respect  to  x  along  contours  of  constant  stream- 
function  (streamlines)  gives 


p  =  F(>V,T)-u^  ^ 

J -00 


dx 


(16) 


where  F  ('1/,T)  is  a  constant  of  integration.  The  density  profile  far  upstream  as  a;  — — cx)  is 
P  ^~p{y)  =P  (^)-  Therefore  F  (T,  T)  =p  (T)  and 


p(x,>V,T)=pm-u^  ^ 

J -00 


dx' 


(17) 


This  implies  that  p  {x,  y,T)  =  p  {x,  T)  and  pt  =  p^i'^T  =  -|-  O  (z/^)  Finally,  combining 

everything  gives 


p  =  p[m)-u  p^  (^) 


'-00 


dx'  +  C>  (u^) 


(18) 


1.1.3  Integration  of  In-Plane  Momentnm  Equation 

We  next  integrate  the  second  equation  in  (13).  Taking  the  y— derivative  of  the  left-hand-side 


gives 


PJ  T)  +  P^^yJ  (T^,  ^ )  -  -  ^^  |  ^  -  =  -  Z/"  (p^^t) 


In  order  to  simplify  notation,  we  dehne 


'y  Ro^  'y 


It, 2  ,  y' 


S  =  p'^yy+P-i,  -4'^  + 


Then 

J(,S,T) 


(19) 


(20) 


=  pJ  {^yy,  ^)  +  Pvl^-/  0^?  +  I’  ^  +  dxU  ^y^yy  +  Px^U  +  I 


p 


p, 


(21) 
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and  the  momentum  equation  can  be  written  as 


J{S,^)  =  \p. 


2^2  ^  y 

2  ^  (3 


X  1^ 


Z/  (p'^yT) 


y  Ro"^ 


(pw) 


(22) 


We  next  make  the  Boussinesq  approximation,  which  can  be  written  formally  as  the  limit  in 
which  the  Boussinesq  parameter,  (3,  approaches  zero.  Physically  this  means  that  the  vertical 
displacement  of  fluid  particles  is  small  compared  with  the  vertical  length  scale  over  which 
the  density  changes  are  signihcant.  In  the  Boussinesq  approximation,  the  vertical  variation 
of  density  scales  as  the  Boussinesq  parameter 

Py=  p^^y  =  Oi/3)  (23) 

and  the  horizontal  variation  of  density  along  a  streamline  can  be  written  using  (15)  as 


I  2PT  2P'^^T 

Pa;  Lt,  =  —u  =  —V  - - - 

\T/  \T/ 


(^2")  =  y‘^l3p  (d')  (d')  +  O  (z/")  (24) 


Then  the  momentum  equation  becomes 


J  (S,  *)  =  N'-‘  (4/)  -  ,,^p(<l)  >1.„T  -  -J^pW 


(25) 


Dehning 


the  momentum  equation  can  be  written  as 


J  (S,  >1.)  = 


(26) 


As  in  the  integration  of  the  density  equation,  we  again  make  the  transformation  {x,y,T)  — > 
{x,  4/,T)  and  assume  the  streamlines  are  not  vertical  throughout  the  domain  (T^  ^  0)Then 
the  momentum  equation  can  be  written 


Cl  -  2p  PW 
P  R  Rq^ 


and  integrated  along  streamlines  to  give 

S  =  Gi^,T)-P  r  ^ 

J -oo 


dx'  — 


Rd^ 


'  -OO  ^y 


dx' 


(27) 


(28) 


The  constant  of  integration,  G  (T,  T),  is  determined  by  the  behavior  of  the  flow  far  upstream 
as  a;  — >•  —oo.  Far  upstream. 


S  —p  (T)  (T)  T  as  a;  — —oo 


G'(T)  =  -p(T)iV2(T)T 


and  the  solution  to  the  momentum  equation  is 


J -OQ  ^ y 


dx  — 


Ro^ 


'  -oo  ^  y 


dx' 


(29) 


Substituting  the  dehnition  of  S  given  in  (20)  and  utilizing  the  Boussinesq  approximation 
gives 


%y  +  (VI/)  (d/  -  //)  = 


Ro^ 


/ 

1-0. 


dx'  +  u^H  +  0  (z/^) 


(30) 


where 
H  = 


d 

d'^ 


^yT^dx'  +  N^{^) 


d 

d'^ 


r 

/ 

'-00  ^y 


dx'  -  yN^  (^) 


d 

(Od' 


r  ^ 

1-00  '^y 


dx' 


We  rewrite  the  hrst  term  on  the  right-hand-side  of  (30)  to  give 

1  d  , 


d/,,  +  (d/)  (d/  -  //)  = 


/?o2 


dx'  +  u^H  +  O 


(31) 


1.1.4  Integration  of  Transverse  (Spanwise)  Momentum  Equation 

We  hnally  integrate  the  equation  governing  momentum  in  the  direction,  given  by  the 
fourth  equation  in  (13).  Using  the  same  coordinate  transformation  as  above,  again  assuming 
d/  ^  0  everywhere,  we  rewrite  this  equation  as 

Wxlvjy  d'y  =  d'y  -  1  -  u^wt  (32) 


and  integrate  to  obtain 

dx  —  f  ^  dx'  (33) 

^  J  — OO  ^  y  ijf 

where  the  constant  of  integration  has  been  set  to  zero  because  tc— vOasa;— >■— cx).  The  first 
term  can  be  rewritten  to  give 


d 


(9^ 


d/ 

-oo  ^  y 


dx' 


(34) 


For  the  remainder  of  the  report,  we  shall  focus  on  the  results  in  the  limit  of  steady  flow  over 
topography.  The  governing  equations  for  steady  flow  are 
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subject  to  the  boundary  condition  along  the  topography 


'h  =  0  on  y  =  eh  (x)  (36) 

and  the  radiation  condition  as  ?/  — >  cx)  ensuring  energy  propagates  outward  away  from  the 
topography. 


1.2  Linear  Solution  for  Uniformly  Stratified  Flow 

In  the  limit  when  the  topography  amplitude  is  small,  e  1,  the  governing  equations  become 
linear  because  the  streamline  displacement,  77  =  y  —  'h  ~  e  scales  with  e.  Thus  derivatives 
with  respect  to  the  streamfunction  can  be  rewritten  as  derivatives  with  respect  to  y  and  the 
linearized  equations  become 


T,,  +  (y)  (4<  -  y) 


1  d 
Ro^  dy 


w  = 


A 

dy 


(T  —  y)  dx' 


(37) 


Dehning  the  streamfunction  pertnrbation,  7/’,  as 


%l)  =  ^  —  y  (38) 

and  combining  equations  gives 

1  52  rx  nx' 

^yy  +  (y)  ^  =  ~^Qy2  J  J  i^dx'dx'  (39) 

Taking  two  derivatives  with  respect  to  x  gives 

(ijyy  +  {y)  ij)  =  0  (40) 

snbject  to  the  linearized  boundary  condition 


tjj  =  —eh  {x)  on  7/  =  0 

Dehning  the  Fourier  transform  in  the  horizontal  and  its  inverse  as 

^  poo 


f{k)=R{f{x)}  = 


271 


f  {x)  e-^^^dx 
f  {k)  R^^dk 


f  (x)  ^  {!  [k)}  ^  I 

J  —00 

the  transformed  differential  equation  becomes 

(k‘^  -  -dyy  +  (y)  ^  =  0 


(41) 


(42) 


(43) 
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subject  to  the  boundary  condition 


-ifj  =  —eh  {k)  on  y  =  0 

and  the  radiation  condition  as  y  — cx),  where  j3  =  ^.  The  general  solution  is 


(44) 


Ip  =  23?  <  /  a{k)  exp  ( i 


k 


y  e  dx  +  /  h{k) exp  -i 


k 


y/k‘^  - 


y  e^^^dx)  (45) 


In  general,  the  radiation  condition  requires  disturbances  to  radiate  outward  from  their  source 
as  time  increases.  Here  the  disturbances  must  radiate  upward  away  from  the  topography. 
The  radiation  condition  can  be  formulated  by  requiring  the  vertical  component  of  group 
velocity  to  be  positive.  The  linear  dispersion  relation  is 


iXJ  = 


and  the  vertical  component  of  the  group  velocity  is 

,  doj  sgn  (km)  k^ 


CgL.  — 


^  dm  m'^  y/W^r~lPm? 


(46) 


(47) 


Therefore  the  radiation  condition  requires  the  signs  of  the  vertical  and  horizontal  wavenum¬ 
bers  to  be  the  same. 

>  0  ^  sgn  {km)  >  0 

In  addition,  the  solution  must  be  bounded  as  cx),  which  requires  h{k)  =  0  when 

\k\  <  13.  Applying  the  radiation  condition  and  boundedness  in  conjunction  with  the  boundary 
condition  at  y  =  0  given  by  (44)  yields  the  solution 


tp  (x,  y)  =  — 2e3?  ■ 


exp 


-ky 


h  {k)  e^'^^dk  + 


exp 


iky 


y/W^ 


h  {k)  d^^dk 


(48) 


The  equation  governing  the  spanwise  velocity  in  the  linear  limit  is  Wx  =  pJy  ci-nd  in  Fourier 

space  w  =  Using  the  solution  for  the  streamfunction  given  in  (48)  gives  the  solution 

for  the  spanwise  velocity  as 


w  {x,  y)  =  — 2e3?  ■ 


ih{k)~ 


Akx 


exp 


-ky 


y/f3‘^  -  k‘^  ^  W(3‘^-k'^ 


dk  -\- 


h  (k) 


:  exp 


iky 


y/W^  ^  \y/k‘^-f3‘^ 


dk 


(49) 

In  both  (48)  and  (49)  the  integral  over  wavenumbers  0  <  k  <  (3  is  evanescent  in  the  vertical 
and  is  the  contribution  dominated  by  rotation.  The  integrals  for  wavenumber  k  >  j3  are 
oscillatory  in  the  vertical.  These  are  the  relatively  short  wavelengths  for  which  rotation  is 
negligible. 
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Figure  3:  Streamlines  for  linear  solution  of  stratified  flow  over  topography  in  the  presence  of  rotation  for 
Ro  =  5.  The  values  of  the  Coriolis  parameter,  Brunt-Vaisala  frequency,  freestream  velocity,  mountain  width 
and  mountain  height  for  this  simulation  were  /  =  1.0  x  10  ^  sec  iVo  =  1.0  X  10  ^  sec  C/q  =  10  m/s, 
Lq  =  20  km  and  ho  =  0.75  km,  respectively.  The  topography  shape  is  Gaussian  (9). 

Figures  3  and  4  show  representative  solntions  for  the  linear  streamlines  and  spanwise  velocity, 
respectively.  The  solutions  shown  are  for  a  Rossby  number  of  5  (/9  =  |).  The  streamlines 
near  a;  =  0  above  the  topography  are  dominated  by  stratihcation  effects,  whereas  further 
downstream  of  the  topography  the  effect  of  rotation  is  evident  in  the  the  longer  wavelength 
inertial-gravity  waves.  The  spanwise  velocity  is  a  purely  rotational  effect  which  is  not  present 
in  the  limit  as  Ro  oo.  Note  that  the  hrst  streamline  in  Figure  3  does  not  coincide  with 
the  topography  because  the  linear  boundary  condition  has  been  applied. 

1.3  Inner  Nonlinear  Solution  and  Lack  of  Uniform  Validity 

While  the  above  solution  accounts  for  rotation,  it  neglects  the  effect  of  nonlinearity.  In  order 
to  incorporate  nonlinearity,  we  return  to  the  original  governing  equations  of  (35).  In  the  case 
of  weak  rotation,  the  Rossby  number  is  large  and  right-hand-side  of  the  second  equation  in 
(35)  is  a  small  correction.  As  a  hrst  approximation,  we  follow  Smith  [21]  and  examine  the 
solution  to  the  nonlinear  problem  with  weak  rotation  using  a  naive  perturbation  expansion 
of  the  form 


^  =  ^(0)  +  +  ^4^(2)  ^  . 

(50) 

Then  the  leading  order  equation  is 

-y  =  0 

(51) 

subject  to 

=  0  on  y  =  eh  (x) 

(52) 
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Figure  4:  Spanwise  velocity,  w  [m/s],  for  linear  solution  of  stratified  flow  over  topography  in  the  presence  of 
rotation  as  given  in  (49)  for  Ro  =  5.  The  parameters  are  the  same  as  those  for  Figure  3. 


and  the  radiation  condition  as  y  ^  oo.  This  is  the  classic  Long’s  problem  [14]  in  which 
the  nonlinear  governing  eqnations  are  rednced  to  a  linear  eqnation  with  nonlinear  bonndary 
condition.  The  solntion  can  be  written  by  decomposing  the  streamfunction  as  above  into  its 
nndistnrbed  component,  y,  and  a  streamfunction  disturbance,  ip. 

+  (53) 

Then  the  solution  becomes 

.0(0)  —  2  (a  (x)  cos y  —  b{x)  sin y)  (54) 

The  radiation  condition  requires 

b  (k)  = —i  sgn  ka  (k)  b  =  —H{a}  (55) 

where  Ti  {a}  is  the  Hilbert  transform  of  a  (x).  Finally,  the  boundary  condition  requires 

a  {x)  cos  (eh)  +  Ti.  {a}  sin  (eh)  =  (x)  (56) 

which  must  be  solved  for  a  (x). 


The  leading  order  equation  for  the  spanwise  velocity  is 

d 


w 


(0)  =  C  (T)  + 


clvl>(0) 


(57) 
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where  the  integral  is  evaluated  along  the  contours  of  constant  Substituting  the  leading 
order  solution  for  gives 


y;(0) 


C'(^)  +2 


d 

,9^(0) 


a  (x')  cosy  —  b  {x')  sin|/]|^(o)  dx' 


(58) 


We  wish  to  examine  the  asymptotic  behavior  of  the  solution  for  the  streamfunction  and 
spanwise  velocity  in  the  limit  as  x  ^  ±cx).  In  order  for  the  Fourier  transform  of  the  topog¬ 
raphy  to  converge  in  the  linear  solution,  it  is  necessary  that  h{x)  =  o  as  x  ^  ±cx).  In 
this  limit  when  x  is  large,  the  amplitude  is  small  and  the  solution  is  linear.  Thus  we  have 


a  ~ 


b  =  —Ti  {a}  =  —i 


2  1  r°° 

d  (k)  e^'^^dk  +  c.c.  ~  -a  (0)  =  —  /  a  (x')  dx' 
X  nx 


>  as  a; 


±CX) 


(59) 

Since  a{x)  b  {x)  as  a;  — >  ±cx),  the  behavior  of  the  leading-order  streamfunction  solution 
as  a;  — >  ±cx)  is 


\]>(o)  ^  y - Q  ag  2;  j-Qo  (60) 

X 

which  decays  to  the  undisturbed  solution,  T  =  y,  as  a;  — >  ±cx).  In  the  limit  as  a;  — ±cx),  the 
spanwise  solution  behaves  as 


Q  roo  _ 2 

~  2—  /  — a  (0)  sinyda;' -I- C  ~  — 4a  (0)  In  bl  cosu  as  a;  — ±cx) 

dy  Jo  X  V  y 

Thus  the  spanwise  velocity  diverges  logarithmically  far  from  the  topography,  and  the  naive 
perturbation  solution  of  (50)  is  not  uniformly  valid,  which  agrees  with  the  hndings  of  Smith 
[21].  As  Smith  suggests,  there  are  really  two  horizontal  length  scales  present.  The  first  is  the 
width  of  the  mountain  used  above,  and  the  second  is  the  length  over  which  rotational  effects 
act,  and  is  Lrot  ~  f-  The  procedure  of  matched  asymptotic  expansions  can  be  used  to 
account  for  these  two  distinct  length  scales  and  the  logarithmic  divergence  of  the  spanwise 
velocity.  In  this  matching,  the  domain  is  divided  into  an  inner  solution  for  a;  — 0  and 
an  outer  solution  as  |a;|  — cx).  The  above  naive  perturbation  solution  represents  the  first 
approximation  for  the  inner  solution,  for  which  rotation  is  a  weak  perturbation.  Far  from 
the  topography  as  length  scales  become  large  and  Coriolis  acceleration  has  sufficient  time 
to  act,  rotation  is  of  leading  order  importance.  However,  assuming  the  topography  decays 
as  h  =  the  disturbances  far  from  a;  =  0  are  small  and  the  solution  becomes  linear, 

allowing  the  application  of  the  linear  boundary  condition  and  superposition.  Figure  (5)  shows 
the  effect  of  rotation  in  the  linear  limit.  Near  the  topography,  the  solutions  with  and  without 
rotation  are  nearly  identical.  Rotational  effects  become  important  further  downstream.  This 
provides  a  justihcation  for  the  approach  of  matched  asymptotic  expansion,  wherein  rotation 
is  weak  in  the  inner  solution  and  is  of  leading  order  importance  in  the  outer  solution. 
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Figure  5:  A  comparison  of  streamlines  showing  the  effect  of  rotation  for  the  linear  solution.  Downstream 
of  the  topography  rotational  effects  become  more  important  as  the  Coriolis  force  has  more  time  to  act. 
Far  downstream  the  inertia-gravity  waves  decay  as  energy  is  radiated  upward.  The  values  of  the  Coriolis 
parameter,  Brunt-Vaisala  frequency,  freestream  velocity,  mountain  width  and  mountain  height  for  Ro  =  10 
were  /  =  1.0  x  10“^  sec“^,  =  1.0  x  10“^  sec“^,  I/q  =  10  m/s,  Lq  =  10  km  and  =  0.75  km,  respectively. 
For  the  linear  limit  of  Long’s  solution  the  parameter  values  were  the  same,  with  the  exception  that  the 
Coriolis  parameter  was  set  to  zero.  The  topography  shape  is  Gaussian  (9). 
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1.4  Outer  Solution 


Far  from  the  center  of  the  topography,  the  solution  is  linear  and  the  governing  equations  are 
given  by  (40)  and  (1.2)  subject  to  the  boundary  condition  given  in  (41).  In  order  to  account 
for  the  leading-order  importance  of  rotation,  the  horizontal  coordinate  and  wavenumber 
should  be  scaled  by  f3  as 

X  =  (3x,  k  =  (62) 

P 

and  the  streamfunction  and  velocity  variables  as 

-ip  =  u  =  1  +  f3u,  V  =  w  =  w  (63) 

where  u  =  ipy  and  v  =  —'ipx-  Then  the  governing  equations  become 


+  N  Pj]  _  +  Pjyy  =  Q 


Wx  =  'Py 

Far  from  the  topography  in  the  stretched  coordinate,  x,  the  topography  appears  to  be  a 
Dirac  distribution.  Then  the  linearized  boundary  condition  along  y  =  0  in  this  outer 
region  can  be  written  as 

pj  =  pp)  =  —2'xeh  (0)  5  ix)  on  y  =  0  (65) 

where  the  factor  27reh  (0)  is  the  area  under  the  topography  curve.  Then  the  solution  to  these 
is  found  via  Fourier  transforms  to  be 

Ip  =  —2eh  (0)  (Ji  -|-  I2) 


w  =  -2eh  (0)  (Is  -h  h) 


where 


Ji  =  3? 


'l-k^ 
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In  order  to  match  this  solution  with  the  inner  solution  near  the  topography,  it  is  necessary 
to  examine  (66)  in  the  limit  as  x  — >  0.  The  limit  of  this  outer  solution  as  x  — >  0  is 


=  27re/i  (0)  siny 

The  details  of  the  evaluation  of  these  limits  is  found  in  Appendix  A. 


1.5  Matching 

In  order  to  match  the  inner  and  outer  solutions,  one  must  return  to  the  inner  solution 
and  evaluate  its  limit  as  x  ^  oo  more  precisely.  The  inner  solution  for  the  streamfunction 
disturbance  was  given  in  (54).  Since  the  integrals  involved  in  evaluating  w  are  along  contours 
of  constant  streamfunction,  it  is  useful  to  rewrite  the  solution  in  terms  of  harmonic  functions 
of  streamfunction,  T,  as 

OO 

•^(0)  =  —  y  =  2  a  (x)  cos  —  6  (x)  sin  ^  sinnT^°^  +  cosn4/*-°^  (69) 

n=l 

where  and  Sn  are  coefficients  arising  from  nonlinearity  and  can  be  determined  analytically 
by  expanding 

OO 

=  2  (a  cos y  —  bsiny)  =  2  a  (x)  cos  —  b  (x)  sin  ^  sin +  Sn  cos 

n=l 

The  hrst  few  terms  can  be  shown  by  hand  to  be 

ri  =  b  (^a'^  +  b^^  ,  Si  =  —a  +  b^^ 

r2  =  2  {a?  —  b^^  ,  S2  =  -^ab 

rs  =  36  —  b^^  ,  S3  =  3a  (sb^  — 
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Substituting  the  above  expansion  into  the  leading  order  equation  for  the  spanwise  velocity 
given  in  (57),  the  spanwise  velocity  is 


sin\h*^°^  /  a  (x')  dx' +  cos  /  b{x')dx' 


E 

n=l 


n 


cosn\l'*'°^  /  (x')  dx' —  sinn\h‘'°^  /  Sn{x)dx 


Then  the  jump  in  spanwise  velocity  across  the  inner  solution  is  found  to  be 


w 


(0) 


n=l 


Details  of  this  computation  are  found  in  Appendix  B. 

In  summary,  we  have  thus  far  obtained  an  inner  solution  of 

CXD 

^(0)  =  2  cos  —  b  sin  ^  (jn  sin  +  Sn  cos 


n=l 


W 


(0)  ^_2 


sin\h‘'°^  /  a{x')  dx' +  /  b{x)dx 


E 

n=l 


n 


cosnT^°^  /  (x')  dx' —  sinn\h*^°^  /  Sn{x')dx' 


+ 


c 


with  asymptotic  behavior 

^(0) 


X 


■u;(o)  ^  — 4a  (0)  cos|/ln  |a;|  +  C 


\x\  oo 


w 


(0) 


oo 

—  OO 


=  — 47rd  (0)  sin  y  +  27r  ^  n  [fn  (0)  cos  ny  —  Sn  (0)  sin  ny] 


n=l 


and  an  outer  solution  of 


^  = -2eh(0)3fJ{/i  +  /2} 

w  =  — 2e/i  (0)  3?  {/a  +  J4} 


(70) 


00 

— 47rd  (0)  sin  +  27r  ^  n  (jn  (0)  cos  (0)  sin  (71) 
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with  asymptotic  behavior 


Ip  ~  2eh  (0) 


sin|/ 

X 


^  ^  0 


TT 


w  ~  2e/i  (0)  cos?/ln  |5;|  —  —  sgn  xsiny 


wIq-  =  27reh  (0)  siny 

At  present,  the  inner  and  outer  solutions  do  not  match.  Thus  it  is  necessary  to  superimpose 
a  second  linear  solution  in  the  outer  domain  upon  the  current  outer  solution  in  order  to 
satisfy  the  matching  conditions  as  5;  — 0.  The  derivation  is  found  in  Appendix  C  and  the 
result  is 

OO  “I  oo 

4a(0)-5]4(0)  3?{/i  +  /2}  +  E4"^ 


w  = 


Aa  (0)  -Y.Sn  (0) 

n=2 

oo 

Ad  (0)  -  E  (0) 


n=2 


(72) 


n=2 


n=2 


where 


=  —2ti  cos  knX  [vn  (0)  sin  ny  +  (0)  cos  ny]  H  [x)  +  2s„  (0)  ,  n  >  1  (73) 

where 


^  =  A  ■ 


j(n)  —  cj: 

Jg  —  ^ 


ke 


ikx 


U2  _  U2 

tv  tv^ 


exp 


ky 


oo  ^glfcX  /  j^y 

exp  i- 


k^  -  kl 


dk  \  ,  n  >  1 


dk  \  ,  n  >  1 


(74) 


This  results  contains  the  proper  jump  condition  for  higher  harmonics,  n  >  2.  However,  the 
hrst  harmonic  in  the  outer  domain  still  does  not  match  the  hrst  harmonic  in  the  inner.  In 
order  to  correct  that,  we  add  another  linear  solution  to  the  solution  in  the  outer  domain, 
such  that  the  the  hrst  harmonic  also  satishes  the  jump  condition.  The  hnal  result  is 


tp  = 


w  = 


Ad  (0) 

n=2 

oo 

Ad  (0)-J2Sn  (0) 


3?  {/i  +  /2}  +  E  4"^  +  2  (0)  (0) 


n=2 


71) 


(75) 


n=2 


^  {/3  +  /4}  +  E  +  2  fsi  (0)  +  h  (0)  w, 


(1) 


n=2 
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where 


=h  +  l7 


1.6  Analytic  Results 

Figures  6  and  7  show  the  streamlines  for  the  linear  and  nonlinear  cases,  respectively  for  a 
Gaussian  mountain  (9)  of  amplitude  e  =  0.75  and  a  Rossby  number  Ro  =  3.  The  solution 
is  only  shown  downstream  of  the  topography  where  the  outer,  nonlinear  analytic  theory  is 
valid  for  a;  3>  0.  Qualitatively  the  linear  and  nonlinear  streamlines  are  very  similar.  This  is 
consistent  with  the  numerical  simulations  of  Triib  and  Davies[22]. 

Figures  8  and  9  are  spanwise  velocity  analogs  of  Figures  6  and  7,  respectively,  showing 
the  spanwise  velocity  far  downstream  of  the  topography.  Here  the  linear  and  nonlinear 
solutions  are  qualitatively  quite  different.  While  the  linear  spanwise  velocity  has  relatively 
long  contours  of  constant  velocity  running  primarily  from  a  high  altitude  upstream  to  a  low 
altitude  downstream,  the  nonlinear  spanwise  velocity  exhibits  much  more  structure  as  shown 
by  the  closed  contours  of  spanwise  velocity  in  Figure  9.  This  is  indicative  of  the  nonlinearity 
in  the  inner  solution,  which  forces  higher  harmonics  in  the  outer  solution.  These  higher 
harmonics  propagate  far  downstream  of  the  topography  and  do  not  dissipate  because  the 
analytic  model  does  not  include  viscous  and  turbulent  dissipative  effects. 

1.7  Computational  Approach 

In  addition  to  an  analytic  theory  for  the  examination  of  rotation  effects  on  nonlinear  flow 
over  topography,  a  computer  code  has  been  developed  for  the  fully  numeric  simulation  of 
these  waves.  The  code  solves  Euler’s  equations  in  a  rotating  coordinate  system  as  given 
in  (1)  using  a  projection  method  for  density-stratified  fluid  based  on  the  work  of  Bell  and 
others  [2,  3]  as  implemented  by  Skopovi  [18,  19,  20]. 
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Figure  6:  Streamlines  from  linear  analytic  solution  far  downstream  of  Gaussian  topography  for  e  =  0.75  and 
Ro  =  3.  The  values  of  the  Coriolis  parameter,  Brunt-Vaisala  frequency,  freestream  velocity,  mountain  width 
and  mountain  height  for  this  simulation  were  /  =  1.0  x  10  ^  sec  fVo  =  1.0  X  10  ^  sec  C/q  =  10  m/s, 
Lq  =  33.3  km  and  ho  =  0.75  km,  respectively.  The  topography  shape  is  Gaussian  (9). 
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Figure  7:  Nonlinear  analytic  streamfunction  in  the  outer  domain  given  by  the  matched  asymptotic  solution. 
The  parameters  used  are  the  same  as  those  for  Figure  6. 
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Figure  8:  Spanwise  velocity  from  linear  analytic  solution  far  downstream  of  Gaussian  topography  for  e  =  0.75 
and  Ro  =  3.  The  values  of  the  Coriolis  parameter,  Brunt-Vaisala  frequency,  freestream  velocity,  mountain 
width  and  mountain  height  for  this  simulation  were  /  =  1.0  x  10  sec  ^,fVo  =  1.0x10  ^  sec  ^,[/o  =  10m/s, 
Lq  =  33.3  km  km  and  ho  =  0.75  km,  respectively.  The  topography  shape  is  Gaussian  (9). 
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Figure  9:  Nonlinear  analytic  spanwise  velocity  in  the  outer  domain  given  by  the  matched  asymptotic  solution. 
The  parameters  used  are  the  same  as  those  for  Figure  8. 
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1.7.1  Overview  of  Projection  Method 

One  of  the  difficulties  in  the  modeling  of  incompressible  flow  is  that  the  time  derivative  of 
pressure  does  not  appear  explicitly,  so  it  is  not  straightforward  to  write  a  time-evolution 
equation  for  it  and  integrate  in  time,  as  done  with  velocity  and  density.  One  approach  to 
resolve  this  is  known  as  the  projection  method  wherein  the  velocity  is  evolved  in  time  using 
a  flnite-difference  approximation,  and  the  time  derivative  of  velocity  is  given  by  all  the  terms 
in  the  Navier-Stokes  equation  with  the  exception  of  the  pressure  gradient  term.  For  a  simple 
forward  Euler  time  discretization,  this  could  be  written  as 

u*  =  (77) 

where  u*  is  an  initial  approximation  to  the  velocity  vector  at  the  next  time  step,  and  f  is 
the  vector  of  forces  acting  on  the  fluid  including  all  stresses  with  the  exception  of  pressure. 
This  new  approximation  of  velocity,  u*  will  not  be  divergence-free  in  general,  which  violates 
conservation  of  mass.  In  order  to  make  the  new  velocity  held  divergence- free,  the  projection 
step  is  used,  whereby  u*  is  projected  onto  a  divergence-free  vector  space.  In  time  the  velocity 
is  advanced  as 

u(n+i)  =  _  AtVp  (78) 

Taking  the  divergence  of  this  gives 

V  ■  =  V  ■  u*  - 

The  left-hand-side  of  this  is  zero  in  order  for  the  flow  to  be  incompressible,  and  the  appro¬ 
priate  pressure  is  determined  by  solving  the  Poisson  equation 

V'p  =  ■  u*  (79) 

Once  the  pressure  is  known,  its  gradient  can  be  determined  and  the  velocity  can  be  updated 
using  (78). 

The  Euler  equations  are  inherently  nonlinear  because  of  the  advective  terms.  In  this  nu¬ 
merical  model,  the  advective  terms  are  treated  using  a  finite-volume  approach  in  order  to 
conserve  momentum.  The  second-order  Godunov  method  is  used  to  approximate  fluxes 
across  finite-volume  boundaries.  The  implementation  follows  the  work  of  Bell  and  others 
[2,  3]  and  is  the  approach  used  by  Skopovi  [18,  19]. 


1.7.2  Viscous  Boundary  Layers 

Throughout  the  analytical  derivation,  the  radiation  condition  was  utilized  to  ensure  energy 
propagates  upwards  and  downstream,  away  from  the  topography.  The  radiation  condition 
was  required  for  the  steady-state  solution.  If  an  initial  value  problem  is  solved,  as  is  the 
case  for  these  fully  numeric  simulations,  a  radiation  condition  is  not  needed  because  energy 
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will  naturally  propagate  outward  from  the  topography.  However,  another  difficulty  arises 
because  the  computational  procedure  is  inherently  limited  to  a  hnite  domain.  Thus  the 
issue  arises  of  what  to  do  when  the  wave  energy  reaches  the  boundary  of  the  computational 
domain.  Methods  of  canceling  the  wave  energy  have  been  proposed  for  a  monochromatic 
spectrum  [7].  However,  in  the  current  problem,  the  spectrum  is  not  monochromatic.  Skopovi 
[18,  19,  20]  has  used  viscous  boundary  layers  to  dissipate  the  wave  energy  before  it  reaches 
the  boundary  of  the  domain. 


In  this  approach,  the  Euler  equations  of  (1)  are  modihed  to  include  a  Newtonian  viscous 
dissipation  term  as 


/Du 

[Im 


-Px  +  p 


d‘^u 


Dv 

=  -Py-  gp  +  P 


'd^v  d^v\ 
dx^  dy"^) 


P 


-Pz  +  P 


(80) 


V  •  u  =  0 


Rp 

Dt 


=  0 


where  p  =  p{x,  y)  is  a  nonuniform  coefficient  of  viscosity.  Note  that  this  formulation  does 
not  agree  with  physical  formulation  of  flow  in  a  held  of  nonuniform  viscosity,  where  the 
gradient  of  the  stress  tensor  appears  on  the  right-hand-side  and 


V  •  (/iVu)  7^  /iV^u 


Nevertheless,  the  viscous  boundary  layers  are  inherently  non-physical  and  our  only  concern 
is  the  prevention  of  rehections  at  the  boundary.  Thus  the  choice  of  viscous  layer  thickness 
and  viscosity  prohle  is  empirical.  Following  the  work  of  Skopovi,  we  choose  a  sinusoidal 
prohle  for  the  viscosity  in  the  viscous  layers.  In  the  viscous  layer  at  the  top  of  the  domain, 
the  viscosity  is  given  as 


.  ^  y-yo\  ^  ^  ^ 

p  =  po  sm  - -  ,  yo<y  <Y 

\2Y  -yoj 

where  y  is  the  vertical  coordinate,  yo  the  height  at  at  which  the  viscous  boundary  layer 
begins,  and  Y  the  height  of  the  top  of  the  computational  domain.  An  analogous  formula 
is  implemented  for  the  viscous  boundary  layers  at  the  left  and  right  edges  of  the  computa¬ 
tional  domain,  and  the  reader  is  referred  to  Skopovi  [19]  for  more  detail.  A  representative 
distribution  of  viscosity  is  shown  in  Figure  10  where  x  and  y  are  nondimensional  coordinates 
used  in  the  computations. 
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Figure  10:  Contours  of  viscosity  for  viscous  boundary  layers  used  to  damp  and  absorb  waves  at  boundaries 
of  domain.  The  viscous  boundary  layers  are  used  as  an  artificial  radiation  condition. 
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The  inclusion  of  viscous  terms  in  an  explicit  temporal  discretization  results  in  a  more  severe 
stability  constraint,  requiring  At  ~  (9  {Ax‘^,Ay‘^)  as  opposed  to  the  stability  constraint  for 
inviscid  flow  of  At  ~  O  {Ax,  Ay).  In  order  to  eliminate  the  need  for  an  extremely  small 
timestep,  the  viscous  terms  are  treated  implicitly  by  solving 

u*  =  +  At/iV^u* 


using  Cholesky  factorization. 

1.8  Computational  Results 

Figures  11  and  12  show  the  streamlines  and  spanwise  velocity,  respectively,  evaluated  using 
the  fully  numeric  simulation  for  stratihed  flow  over  topography  with  weak  rotation.  The 
results  are  for  a  nondimensional  mountain  amplitude  of  e  =  0.40  and  the  results  are  scaled 
up  linearly  to  an  amplitude  of  e  =  0.75  for  comparison  with  the  nonlinear  results  shown  in 
Figures  13  and  14.  The  results  are  qualitatively  quite  similar  to  the  linear,  analytic  solution 
shown  in  Figures  3  and  4  with  different  scaling  because  a  Rossby  number  of  5  had  been  used 
for  the  linear,  analytic  results  previously.  The  waves  propagate  upward  and  downstream 
(rightward)  from  the  topography.  The  plotted  solution  has  not  yet  reached  steady-state  far 
aloft  and  downwind  of  the  topography.  The  results  are  plotted  for  this  time  despite  having 
not  achieved  steady-state  because  as  time  increases,  the  viscous  layers  become  saturated  and 
are  no  longer  able  to  dissipate  all  waves  at  the  boundary.  Thus  for  large  times,  spurious 
reflections  from  the  viscous  layers  and  boundaries  are  found  in  the  interior  of  the  domain. 
The  results  shown  in  Figures  11  through  14  have  negligible  reflections  but  have  not  reached 
steady-state  throughout  the  domain. 

Figures  13  and  14  show  the  numerically  calculated  streamfunction  and  spanwise  velocity  re¬ 
spectively.  These  plots  are  analogous  to  those  in  Figures  11  and  12  with  the  only  difference 
being  that  the  topography  amplitude  is  e  =  0.75  for  Figures  13  and  14,  whereas  it  is  e  =  0.40 
for  Figures  11  and  12. 

Qualitatively,  the  numeric  results  with  e  =  0.40  differ  very  little  from  those  with  e  =  0.75. 
While  the  analytic  theory  showed  much  more  structure  in  the  spanwise  velocity  for  the  case 
of  e  =  0.75,  the  fully  numeric  solution  does  not  show  this.  Part  of  the  reason  for  this  is  that 
the  numerical  solutions  have  achieved  steady-state  for  a  very  limited  horizontal  domain.  The 
analytic  results  in  Figure  9  show  closed  contours  of  spanwise  velocity  far  downstream  of  the 
topography,  and  very  few  of  these  contours  are  present  in  the  domain  x  <  900  km,  for  which 
we  have  numeric  results. 


1.9  Discussion 

The  results  of  this  section  show  that,  for  a  single  layer  of  uniformly  stratihed  how,  nonlinear 
interactions  above  the  mountain  generate  higher-harmonic  disturbances  that  propagate  far 
downstream  of  the  topography.  This  is  best  seen  in  comparing  Figure  9  for  nonlinear  span- 
wise  velocity  with  Figure  8  for  the  corresponding  linear  result  of  spanwise  velocity.  While 
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Figure  11:  Streamlines  from  fully  numeric  simulation  for  relatively  weak  nonlinearity.  The  values  of  charac¬ 
teristic  Coriolis  parameter,  Brunt -Vaisala  frequency,  undisturbed  wind  speed,  mountain  width  and  mountain 
height  for  this  simulation  were  /  =  1.0  x  10“"^  sec“^,  Nq  =  1.0  x  10“^  sec“^,  Uq  =  10  m/s,  Lq  =  33.3  km  and 
/iQ  =  0.40  km,  giving  nondimensional  parameters  Ro  =  3,  e  =  0.40  and  fj,  =  ^.  The  streamlines  are  plotted 
at  a  computational  time  of  tcomp  =  200,  corresponding  to  a  nondimensional  time  of  t  =  33.33  or  a  physical 
time  of  tdimensionai  ~  185  hr.  The  amplitude  of  the  streamfunction  disturbance  is  scaled  by  a  factor  of  1.875 
for  comparison  with  results  for  the  larger-amplitude  mountain  below. 
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Figure  12:  Spanwise  velocity,  w  [m/s],  from  fully  numeric  simulation  for  relatively  weak  nonlinearity.  The 
parameters  for  these  results  are  the  same  as  those  for  Figure  11.  The  amplitude  of  the  spanwise  velocity  is 
scaled  by  a  factor  of  1.875  for  comparison  with  results  for  the  larger-amplitude  mountain  below. 
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Figure  13:  Streamlines  from  fully  numeric  simulation  for  mountain  with  nondimensional  amplitude  e  =  0.75. 
The  parameters  are  the  same  as  those  for  Figure  12,  with  the  exception  of  the  topography  height  and 
nondimensional  amplitude,  for  which  the  values  are  ho  =  0.75  km  and  e  =  0.75,  respectively.  Furthermore, 
the  amplitude  of  the  response  is  unsealed. 
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Figure  14:  Spanwise  velocity  from  fully  numeric  simulation  for  relatively  strong  nonlinearity. 
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these  higher  harmonics  are  present  for  a  single  layer  atmosphere,  they  are  relatively  weak 
because  the  nonlinear  interactions  in  Long’s  solution  for  a  single  layer  are  relatively  weak. 
Thus  the  analysis  serves  as  a  proof-of-concept  showing  inertia-gravity  waves  far  downstream 
produced  by  nonlinear  interactions  above  the  mountain,  but  these  waves  have  relatively 
small  amplitude  due  to  weak  nonlinear  interactions.  The  subsequent  section  shows  that  in  a 
two-layer  atmosphere,  the  effect  of  the  tropopause  can  signihcantly  increase  the  amplitude 
of  the  nonlinearly  generated  inertia-gravity  waves  downstream. 

The  numerical  results  in  this  section  are  rather  inconclusive.  The  primary  cause  is  limited 
computational  resources.  The  fully  numeric  simulation  is  constrained  to  a  hnite  domain  and 
the  simulations  are  started  from  a  uniform  velocity  held  and  allowed  to  evolve  in  time.  The 
spatial  and  temporal  domains  shown  in  Figures  11  through  14  are  insufficient  to  conclusively 
resolve  nonlinearities  far  downstream  of  the  topography  at  steady-state.  For  example,  the 
horizontal  domain  for  the  numerical  results  is  limited  to  a;  <  900  km  in  the  numerical  results 
of  Figures  11  through  14.  This  implies  that  the  analytical  theory  provides  a  predictive  model 
yielding  information  unable  to  be  obtained  from  these  numeric  simulations. 

In  summary,  the  results  of  this  analysis  show  that  nonlinearity  in  the  inner  solution  slightly 
affects  the  spanwise  velocity  in  the  outer  solution.  However,  this  forcing  is  relatively  weak 
because  the  nonlinearities  above  the  topography  in  Long’s  solution  for  a  single,  uniformly 
stratihed  layer  are  relatively  weak,  and  the  nonlinear  solution  does  not  differ  appreciably 
from  Long’s  solution  in  the  linear  limit.  Next  we  examine  the  effect  of  the  tropopause  in  a 
model  atmosphere  with  two  layers  of  uniformly  stratihed  how  and  seek  conditions  for  which 
the  nonlinear  forcing  of  the  outer  solution  is  stronger. 
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2  Two  Layers  of  Uniform  Stratification  —  Effect  of  the  Tropopause 
2.1  Governing  Equations 

In  general,  there  is  a  change  in  the  Brunt-Vaisala  frequency  across  the  tropopause,  with 
the  Brunt-Vaisala  frequency  typically  being  smaller  in  the  troposphere  and  larger  in  the 
stratosphere.  While  the  Brunt-Vaisala  frequency  is  not  spatially  uniform  in  either  layer  of 
the  atmosphere,  we  formulate  a  model  to  examine  the  effect  this  change  in  Brunt-Vaisala 
frequency  across  the  tropopause  has  in  conjunction  with  nonlinearity  and  rotation.  The 
effect  of  the  tropopause  for  cases  excluding  rotation  has  been  previously  studied  by  several 
investigators  including  Durran[6]  and  Davis[5].  In  this  investigation,  we  utilize  the  frame¬ 
work  of  matched  asymptotics  developed  above  to  analyze  a  system  with  two  layers  of  uniform 
Brunt-Vaisala  frequency  in  the  presence  of  nonlinearity  and  rotation.  The  model  under  in¬ 
vestigation  is  sketched  in  Figure  15. 


Fo 


Si 


V’2  (x,  y) 


Figure  15:  Physical  system  of  stratified  flow  over  topography  in  the  presence  of  rotation  incorporating  the 
effect  of  the  tropopause. 


The  governing  equations  for  this  system  are  those  derived  in  §1.1  and  given  in  (35)  and 
(36)  where  we  set  the  Brunt-Vaisala  frequency  to  Ni  in  the  troposphere  and  W  in  the 
stratosphere,  with  a  step-change  at  the  tropopause. 
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where  is  the  streamline  on  the  tropopause  separating  the  troposphere  and  stratosphere. 
As  in  the  case  of  a  single  layer  of  stratihed  flow,  the  kinematic  boundary  condition  at  the 
surface  of  the  topography  is  imposed,  implying  that  there  is  no  flow  normal  to  the  surface 
and  resulting  in  the  fact  that  the  topography  is  itself  a  streamline. 

'h  =  0,  y  =  eh  {x)  (82) 

The  radiation  condition  as  y  ^  oo  must  also  be  imposed.  Finally,  we  must  impose  some 
interfacial  conditions  at  the  tropopause.  The  first  is  the  kinematic  condition  requiring  the 
flow  in  the  vertical  direction  to  be  equal  on  either  side  of  the  tropopause,  resulting  in  the 
fact  that  tropopause  itself  is  a  streamline  and  the  streamfunction  on  either  side  must  be 
equal  to  zero. 

'h(i)  =  T(2)  =  T*,  T  =  T*  (83) 

Since  the  mass  flow  rate  between  two  streamlines  is  proportional  to  the  difference  in  stream- 
function,  this  is  consistent  with  the  fact  that  the  troposphere  and  stratosphere  are  in  contact 
with  each  other  and  there  is  no  flow  between  the  two  layers.  This  assumes  the  tropopause 
has  inhnitesimal  thickness,  consistent  with  the  model  proposed  in  (81). 


The  second  interface  condition  is  a  dynamic  condition  coming  from  a  balance  of  pressure  on 
either  side  of  the  tropopause  streamline,  T*. 

Pi  =  P2,  ^ 


which  implies  ui  =  M2  on  T  =  T*  in  the  hydrostatic  limit.  Rewriting  the  horizontal  velocity 
in  terms  of  the  streamfunction  gives  the  second  boundary  condition  at  the  tropopause. 


<9^(1)  ^  <9T(2) 
dy  dy 


(84) 


2.2  Linear,  Uniformly  Valid  Solution 

In  the  limit  of  e  — 0,  the  streamfunction  disturbances  are  small  (of  order  e)  and  surfaces  of 
constant  altitude  approach  surfaces  of  constant  streamfunction. 

'^  =  y  +  iIj,  ifj  =  O  (e) 

Then  the  boundary  conditions  at  the  topography  and  tropopause  can  be  applied  along 
contours  of  constant  altitude.  Decomposing  the  streamfunction  into  its  undisturbed  value, 
y  and  a  streamfunction  disturbance,  ip,  the  governing  equations  become 
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subject  to 


^(1)  =  -eh{x) ,  1/  =  0 
^(1)  =  ^(2),  y  =  H 
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dy  dy 

and  the  radiation  condition.  The  solution  is 


\y  =  H 


^(i)=  23?. 


fl3 


—eh  cosh 


Niky 


+  g  sinh 


—eh  cos 


Niky 

Vk'^  - 


+  g  sin 


Niky 
V/9"  -  k'^ 

Niky 


e^^^dk 


e^^^dk),  y<H 


(  -N,ky 


:  N^ky  \  Jkx, 


dexp  /  cexp^«-^  e-d.  .  y>H 


tc(i)  =  23? 


-iNi 


W(2)  =  23? 


7/92  _  ^2 

/•“  -iiVi 

'y 

i  iN2 


7/3^  - 

eh  sinh 


eh  sin 


iVihi/ 

V(3^-k\ 
Niky 


+  5f  cosh 


y/k‘^-p\ 
Niky 


g  cos 


7/32  -p 

iVih|/ 


^y/k^-f3\ 


e^’^^dk),  y<H 


dexp 


^  -N2ky  , 


7/32  _  ^2 

/  1’  * 


(87) 


32 


with  constants  of  integration  c{k),  d  {k)  and  g  {k)  determined  from  the  matching  conditions 
to  be 


c{k)  = 
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The  details  of  this  solution  are  found  in  Appendix  D. 


(88) 


2.3  Inner  Solution 


The  preceding  solution  neglects  the  effect  of  nonlinearity,  which  is  the  primary  concern  of 
this  study.  In  order  to  account  for  nonlinearity,  a  matched  asymptotic  solution  is  derived 
as  in  the  case  of  a  single  layer  of  uniformly  stratihed  flow.  In  the  inner  solution  near  the 
topography,  treated  as  a  weak  perturbation  to  the  two-layer  analog  of  Long’s  solution.  The 
governing  equations  for  the  inner  solution  are  then 


p)  ^  0,  ^ 


dy'^ 


wW 


y)  =0,  T  >  T* 
dx' 


y 


(89) 


where  the  superscript  (0)  denotes  the  leading-order  solution  and  the  subscripts  denote 
whether  the  solution  is  in  the  troposphere,  (1),  or  the  stratosphere,  (2).  Here  the  stream¬ 
line,  T  =  T*  is  the  streamline  of  the  tropopause,  which  is  not  known  a  priori  and  must  be 
determined  as  part  of  the  solution  of  the  equations. 
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The  problem  of  two  layers  of  uniformly  stratified  flow  in  the  absence  of  rotation  is  an  ex¬ 
tension  of  Long’s  classic  solution  [14]  and  has  been  studied  by  Durran[6]  and  Davis [5].  The 
solution  can  be  formulated  generally  as 


■^(1)  =  a  (x)  cos  (Niy)  -|-  b  {x)  sin  {Niy) ,  T  <  T* 

(90) 

^/J^2)  =  c  (x)  COS  {N2y)  +  d  {x)  sin  (iVsy) ,  T  >  T* 

where  we  have  decomposed  the  streamfunction,  T,  into  the  sum  of  its  undisturbed  value, 
y  and  a  streamfunction  disturbance,  Then  the  kinematic  boundary  condition  at  the 
topography,  y  =  h{x),  requires 


a  cos  {Nih)  +  6  sin  {Nih)  =  —h{x) 


(91) 


and  the  radiation  condition  as  y  — >  cx)  gives 


d{x)  =  H{c{x)}  (92) 

resulting  in  wave  energy  propagating  vertically  upwards.  As  in  the  linear  case,  two  conditions 
are  needed  at  the  tropopause  to  match  the  solutions  in  the  troposphere  and  stratosphere. 
We  can  write  the  curve  of  the  troposphere  as 


^*  =  H +  'ip*  {x)  =  H -f]  {x) 


(93) 


where  H  is  the  undisturbed  height  of  the  tropopause  and  rj  is  the  vertical  displacement  of  the 
tropopause,  giving  the  tropopause  location  as  y  =  H  +  rj  (x),  which  is  also  unknown.  The 
hrst  matching  condition  at  the  tropopause  is  the  kinematic  matching  condition  requiring 
the  streamfunction  in  the  troposphere  and  stratosphere  to  be  equal  at  the  tropopause.  This 
gives 

a  cos  [Ni  {H  +  7])]  +b  sin  [Ni  {H  +  y)]  =  c  cos  [W  {H  +  y)]  +  d  sin  [W  {H  +  y)]  (94) 

The  second  interface  matching  condition  is  the  dynamic  bonndary  condition  reqniring  the 
pressure  to  be  equal  on  either  side  of  the  tropopause.  Integrating  Bernoulli’s  equation  from 
far  upstream  and  noting  that  the  density  in  the  troposphere  at  the  tropopause  is  equal  to 
the  density  in  the  stratosphere  at  the  tropopause,  gives 

U(i)  =  u^2),  y  =  H  +  y  =>  1/ =  if -^  77  (95) 

Substituting  the  general  solntion  of  (90)  gives 

Ni  {— asin  [Ni  {H  +  y)]  +  bcos  [W  (if  -|-  y)]}  =  N2  {— csin  [N2  {H  +  y)]  +  dcos  [N2  {H  +  y)]} 

(96) 

Finally,  the  hfth  equation  needed  is  simply  a  specihcation  of  the  tropopause  displacement, 

y. 

y  =  —t/j*  =  —a  {x)  cos  [Ni  (if  +  y)]  —  b  {x)  sin  [Ni  (if  -|-  77)]  (97) 

Thus  (91),  (92),  (94),  (96)  and  (97)  are  the  hve  eqnations  used  to  determine  a{x),  b{x), 
c{x),  d{x)  and  y{x).  These  are  system  of  nonlinear  equations  becanse  the  bonndary  and 
matching  conditions  are  applied  along  contours  of  constant  streamfunction,  as  opposed  to 
contours  of  constant  altitude.  They  can  be  solved  nnmerically  using  an  iterative  procednre, 
snch  as  Newton-Raphson. 
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2.4  Outer  Solution 

The  outer  solution  is  obtained  by  scaling  the  linear,  uniformly  valid  solution  of  (87)  into  the 
outer  variables  as  done  for  the  single-layer  case  above.  Without  loss  of  generality,  we  can  set 
the  nondimensionalized  Brunt-Vaisala  frequency  in  the  stratosphere  equal  to  unity,  W  =  1, 
and  the  nondimensionalized  Brunt-Vaisala  frequency  in  the  tropopause  to  Ni  =  N.  Then 
the  outer  solution  is 

■^(1)  =  —2  (/i  -I-  12) 


'^{2)  —  —2  (/s  +  74) 

u'(i)  =  —2  -f  IqJ 


(98) 


U'(2)  —  —2  (^ij  +  /s) 

where  the  integrals  are  dehned  in  Appendix  E  and  have  asymptotic  behavior  as  7  — 0  of 


^(1) 

^(2) 

h)(i) 


2eh  (0)  - 


N  sin  Ny 


X  cos2  NH  +  sin^  NH 


2eh  (0)  - 


N 


X  m  cos2  NH  +  sin^  NH 
2eh(0)  V 


sin  NH  cos  {y  —  H)  +  N  cos  NH  sin  {y  —  H)] 


cos2  NH  +  sin^  NH 

TT 


{N  In  \x\  cos  Ny 


sgn  X 


N^  cos  NH  sin  N  {y  —  H)  +  sinNH  cos  N  (y  —  H) 


W(2) 


2eh{0)N 


{In  \x\  [V  cos  Videos  {y  —  H)  —  sin  VTfsin  {y  —  H)] 


V2  cos2  NH  +  sin^  NH 

TT  I 

—  sgn  X  [N  cos  NH  sin  (y  —  H)  +  sin  NH  cos  {y  —  H)]  | 

and  a  jump  in  the  spanwise  velocity  across  the  inner  domain  of 
27reh  (0)  N 


>  ,  X 


(99) 


W(i) 


0+ 


0-  m  cos2  NH  +  sin^  NH  L 


N^  cos  NH  sin  N  {y  —  H)  +  sin  N H  cos  N  {y  —  H) 


W(2) 


27ieh  (0)  N 


0+ 


0-  Ar2  cos2  NH  +  sin^  NH 


[N  cos  NH  sin  (y  —  H)  +  sin  NH  cos  {y  —  H)] 


(100) 
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2.5  Matching 

Rewriting  the  inner  solution  as 


1^(1)  ~  ^  (^)  (^) 

1^(2)  ^  P  (^)  (l/  ~  {P  (2^)}  sin  {y  —  H) 


(101) 


the  asymptotic  behavior  as  a;  — >  cx)  is 


p  (x)  ~  — 
X 


n{p  (x)} 


CN  1 


tan  NH  x 
r  (x)  ~  —eh  (x) 
Cl 

s  (x)  ^  — 


sin  NH  X 


>,  X  ^  00 


(102) 


T]  (x)  ~  —p  (x) 

where  the  constant  C  must  be  determined  numerically.  In  the  linear  limit,  the  constant  C 
approaches 

2Neh  (0)  sin  NH 


Clin  — 


Thus 


m  cos2  NH  +  sin^  NH 
C  sin  Ny 


i>m  ~  c 


sin  NH  X 
cos  (p-Tr)  +  ^;^sin  jy-H) 

X 


>,  X  ^  00 


(103) 


(104) 


In  order  to  properly  match  the  jump  in  spanwise  velocity,  it  is  useful  to  rewrite  the  inner, 
nonlinear  solution  in  terms  of  harmonics  of  the  streamfunction  as 

^  W  =  r  cos  N  -  iP)  +  s  sin  iV  -  h) 

00 

l~  {x)  sin  niV  —  p)  +  rn~  (x)  cos  nN 


n=l 


V(2) 


p  COS  —  H^  +H  {p}  sin  —  H^ 

00 

V  (x)  sinn  —  y)  +  (x)  cos  71  —  y) 


(105) 


n=l 
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where  and  are  the  coefficients  (as  a  function  of  x)  of  the  terms  arising  from  nonlinear 
interactions  in  the  inner  solution.  These  terms  will  give  rise  to  jumps  in  spanwise  velocity  at 
higher  vertical  harmonics  and  are  the  forcing  of  the  higher  harmonics  in  the  outer  solution. 
They  can  be  evaluated  by  orthogonal  decomposition  as 


L  (a^)  = 


N 


TT 


N 


[X)  =  — 
71 

^2Tr 


-y)  sinniVT(°)dT(°) 
-  I/)  cosniVT(°)dT(°) 


C  {x)  =  -  I  (t(°)  -  y)  sinn  -H)d  -  h) 


m+  (x)  =  - 

TT 


'»27r 

-  y)  cosn  -  h)  d  (t®  -  h) 


,  n  =  2,3,4,---  (106) 


and  if  and  mf  can  be  found  by  subtracting  the  appropriate  leading  order  function,  p{x), 
'H{p{x)},  r  (x)  or  s  {x),  from  the  orthogonal  decomposition. 


h  (x) 

=  -s  ix)  H - 

TT 

mj“  {x) 

=  — r  (x)  H - 

71 

If  (x)  =  -p  (x)  + 


r»27r 


TT 


mf  (x)  =  —Id  {p  (x)}  + 
Then  the  spanwise  velocity  is 


-  I/)  cosiVT(°)dT(°) 

-  y)  sin  -  h)  d  -  h) 

•27r 

-  y)  cos  (t'”’  -H)d  (ill”’  -  H 


(107) 


w 


(0) 

(1) 


(0) 

^(2)  = 


=  —  iVsiniVT*-'^^  /  r  (x')  dx' +  iVcosiVT*-®^  /  s{x')dx 

Jo  Jo 

oo  f  nx  nx 

+A^  ^  n  jcosniVT*^®^  y  l~{x)dx—smnN'dj‘^^^j  m“(a;')da;'|  +  C(i) 

nx  px 

sin  j  p  {x)  dx  +  cos  J  H{p{x')}dx' 

OO  (  pX  pX  'j 

+  ^  n  cos  n  “  ^)  /  C {x)d  X  —  sin  n  (*(“l  -H)  m7,{x')dx'  +  C(2) 

(108) 
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Using  the  asymptotic  behavior  of  the  functions  p{x),  'H{p{x)},  r  [x)  and  s{x)  given  in 
(102),  the  asymptotic  behavior  of  the  spanwise  velocity  as  a;  — >  cx)  is 


w 


w 


(0) 

(2) 
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(0) 

(1) 
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CN 


cos  Ny  In  |a;| 


tan  NH 


sin  NH 
cos  {y  —  H)  —  sin  {y  —  H) 


)  ,  X 


In  \x\ 


oo 


(109) 


exhibiting  the  same  logarithmic  divergence  and  lack  of  uniform  validity  as  seen  in  the  single 
layer  of  uniformly  stratified  flow  and  conhrming  the  need  for  a  matched  asymptotic  expan¬ 
sion.  In  addition  to  this  logarithmic  divergence,  the  spanwise  velocity  experiences  a  hnite 
jump  across  the  inner  domain  given  by 


w 


=  —N  sin  Ny  /  r  {x')  dx  -|-  N  cos  Ny  /  s  (x')  dx 


+N  X!  ^  nNy  /  {x)dx'  —  sin  nNy  /  {x)dx  \ 

n=l  {  J —oo  J -oo  ) 

/OO  POD 

p  {x')  dx' +  COS  {y  —  H)  /  H{p{x')}dx' 

OO  J  — OO 

oo  (  POO  POO  'j 

+  '^n<cosn{y  —  H)  /  l^{x')dx'  —  sin n{y  —  H)  /  m'^{x')dx' \  +  C(2) 

n=l  I  J-oo  J -oo  j 

(110) 

In  order  to  match  the  inner  solution,  we  write  the  outer  solution  as  a  multiple  of  the  linear 
outer  solution  given  in  (98)  plus  terms  arising  from  nonlinearity  in  the  inner  solution. 


^(1)  -  (Ji  -h  I2)  +  '4’(l),nonlin  +  ^(l),x 
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n=2 


J,(0)  _  9  4  m  I  r  'i  I  1 
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(111) 


^(ij  -  (Is  -f  Je)  *^(iy 


(1) 


dn) 


n=2 


-2A  (I7  +  Is)  + 


OO 

+  ^4. 


n=2 


(n) 

(2),y 


38 


where  ^2)  satisfy 

=  -27riV^  (n^  -  l)  6'  (x)  |/“  (0)  sinnNy  +  m~  (0)  cosnNy'^ 

+  iV^<h(2)  j  +  =  -27r  (n^  -  l)  S'  (x)  |/+  (0)  sinn  (?/  -  H)  +  (0)  cosn{y  -  H 

(112) 

subject  to  the  boundary  condition  at  the  topography,  the  matching  conditions  at  the  tropopause 
and  the  radiation  condition  as  y  — >  cx). 


2.6  Results  and  Discussion  —  Analytic 

In  this  two-layer  model  of  the  atmosphere,  the  height  of  the  tropopause  affects  the  amplitude 
of  the  waves.  For  certain  tropopause  heights,  the  amplitude  of  the  waves  can  be  signihcantly 
enhanced,  which  is  commonly  referred  to  as  ‘tropopause  tuning’.  Likewise,  detuning  ex¬ 
ists  for  other  tropopause  heights  giving  attenuated  wave  amplitudes.  The  phenomenon  of 
tropopause  tuning  has  been  previously  studied  by  Klemp  and  Lilly [11]  and  others  who  found 
that  the  response  in  the  linear  limit  is  tuned  when  the  tropopause  is  located  at  odd  multiples 
of  half  the  vertical  wavelength  in  the  troposphere. 

In  order  to  investigate  the  phenomenon  of  tropopause  tuning  and  the  effect  it  has  on  inertial- 
gravity  waves  far  downstream  of  the  topography,  we  examined  the  response  of  the  inner 
solution  where  rotation  is  negligible.  The  streamline  slope  is  a  measure  of  the  nonlinearity 
of  the  solution,  where  the  linear  solution  has  slope  of  order  e  — 0  and  the  nonlinear  solution 
has  a  hnite  slope.  Figure  16  shows  the  maximum  streamline  slope  in  the  lower  layer  of  the 
inner  solution  as  a  function  of  tropopause  height.  It  is  seen  that,  for  this  range  of  tropopaue 
heights,  there  are  three  distinct  peaks  at  which  the  streamline  slope  is  a  local  maximum, 
where  the  tropopause  height  is  tuned. 

The  nonlinearity  in  the  inner  solution  manifests  itself  in  the  outer  solution  through  the  non¬ 
linear  coefficients  dehned  in  (106)  and  (107).  In  order  to  verify  that  a  tropopause  height 
with  steep  streamlines  will  provide  the  largest  amount  of  nonlinear  forcing  on  the  outer  so¬ 
lution,  these  nonlinear  coefficients,  and  m^,  were  examined  as  a  function  of  tropopause 

/OO 

m2  (x)  dx  as  a  function  of  the  tropopause 

■OO 

height.  The  coefficient  of  the  second  harmonic  was  chosen  because  quadratic  nonlinearities 
were  larger  than  the  higher  harmonics.  Again  the  three  distinct  local  maxima  are  observed 
at  approximately  the  same  tropopause  heights  as  those  in  Figure  16,  conhrming  that  steep 
streamlines  give  rise  to  larger  nonlinearities. 

Based  upon  the  results  of  the  streamline  slope  in  Figure  16  and  the  coefficients  of  nonlinear¬ 
ity,  of  which  Figure  17  is  representative,  two  differening  cases  were  selected  for  comparison. 
A  relatively  tuned  case  with  a  tropopause  height  of  H  =  12.5  km  was  compared  with  a 
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Figure  16: 
0.375. 


Maximum  slope  of  streamlines  in  inner  solution  as  a  function  of  tropopause  height  for  e 
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Figure  17:  A  representative  coefficient,  (0),  of  nonlinearity  below  the  tropopause  giving  rise  to  forcing 
of  the  outer  solution  by  the  nonlinearity  in  the  inner  solution.  The  n  =  2  mode  is  the  quadratic  nonlinearity 
and  has  the  largest  amplitude  for  this  value  of  e  =  =  0.375. 
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Figure  18:  Streamfunction  for  flow  of  two  layers  of  uniformly  stratified  fluid  with  mean  velocity  Uq  =  10  m/s 
over  a  mountain  of  height  h  =  0.375  km  centered  at  x  =  0  km  with  Rossby  number  Ro  =  3.  Tuned  case 
with  tropopause  height  H  =  12.5  km.  This  case  shows  waves  propagating  far  downwind  of  the  topography 
with  very  little  decrease  in  amplitude. 

relatively  detuned  case  with  tropopause  height  H  =  11.0  km.  Figures  18  and  19  show  the 
analytical  solutions  for  streamlines  and  spanwise  velocity  far  downstream  of  the  topography 
for  the  tuned  case  with  tropopause  height  H  =  12.5  km.  This  is  contrasted  with  the  detuned 
case  shown,  for  which  the  analogous  results  are  shown  in  Figures  20  and  21. 

The  tuned  and  detuned  streamlines  in  Figures  18  and  20,  respectively,  are  qualitatively 
quite  different.  In  the  troposphere,  the  disturbance  is  significantly  stronger  for  the  tuned 
case.  The  waves  in  the  troposphere  propagate  far  downstream  (3000  km)  with  no  signihcant 
attenuation.  The  detuned  streamlines  in  the  troposphere  show  relatively  weak  waves  that 
decay  downstream  of  the  mountain  and  are  all  but  nonexistant  3000  km  downstream  of  the 
topography.  In  the  stratosphere,  the  differences  between  the  tuned  and  detuned  cases  are 
not  as  severe.  Both  exhibit  qualitatively  similar  streamline  patterns,  but  in  the  tuned  case 
the  amplitude  is  larger  further  downstream. 

The  spanwise  velocity  prohles  shown  in  Figure  19  for  the  tuned  case  and  Figure  21  for  the 
detuned  case  exhibit  the  same  qualitative  features  as  the  streamlines.  In  the  tropopause,  the 
tuned  response  is  much  stronger  than  the  detuned  response  and  is  nearly  as  strong  3000  km 
downstream  of  the  mountain  as  it  is  500  km  downstream  of  the  mountain.  The  detuned 
spanwise  velocity  in  the  troposphere  is  weaker  than  the  tnned  response,  and  its  amplitude 
decays  signihcantly  as  the  distance  downwind  of  the  mountain  increases.  In  the  stratosphere, 
the  spanwise  velocity  exhibits  the  closed  contours  characteristic  of  nonlinearities;  whereas, 
the  tnned  response  shows  relatively  long,  smooth  contonrs  of  spanwise  velocity.  This  is 
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Figure  19:  Contours  of  transverse  velocity  for  flow  of  two  layers  of  uniformly  stratified  fluid  with  mean 
velocity  C/q  =  10  m/s  over  a  mountain  of  height  h  =  0.375  km  centered  at  x  =  0  km  with  Rossby  number 
Ro  =  3.  Tuned  case  with  tropopause  height  H  =  12.5  km.  This  case  shows  waves  propagating  far  downwind 
of  the  topography  with  very  little  decrease  in  amplitude. 
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Figure  20:  Streamfunction  for  flow  of  two  layers  of  uniformly  stratified  fluid  with  mean  velocity  Uq  =  10  m/s 
over  a  mountain  of  height  h  =  0.375  km  centered  at  x  =  0  km  with  Rossby  number  Ro  =  3.  Detuned  case 
with  tropopause  height  H  =  11  km.  This  response  is  relatively  weak  and  decays  far  downstream  of  the 
topography. 
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Figure  21:  Contours  of  transverse  velocity  for  flow  of  two  layers  of  uniformly  stratified  fluid  with  mean 
velocity  Uq  =  10  m/s  over  a  mountain  of  height  h  =  0.375  km  centered  at  a:  =  0  km  with  Rossby  number 
Ro  =  3.  Detuned  case  with  tropopause  height  H  =  11  km.  This  response  is  relatively  weak  and  decays  far 
downstream  of  the  topography. 

similar  to  the  observations  in  the  single  layer  of  nniformly  stratihed  flow  for  the  linear 
solntion  shown  in  Fignre  8  and  the  nonlinear  solntion  shown  in  Fignre  9. 

2.7  Results  and  Discussion  —  Numeric 

In  addition  to  the  above  analytic  investigation  of  tropopanse  height,  a  companion  nnmerical 
simnlation  was  rnn  to  investigate  the  combined  effect  of  the  tropopanse  and  nonlinearity. 
As  stated  above,  the  slopes  of  the  streamlines  in  the  inner  solntion  along  with  the  nonlin¬ 
earity  coefficients  indicate  that  the  tropopanse  is  relatively  tnned  at  heights  of  if  ~  6.5  km, 
10  km  and  12.5  km.  Becanse  the  inner  solntion  is  periodic,  the  effect  of  tropopanse  height  is 
periodic  with  a  period  eqnal  to  the  vertical  wavelength  in  the  troposphere.  Thns  the  tnning 
seen  at  H  =  12.5  km  investigated  analytically  above  shonld  be  approximately  eqnal  to  the 
tnning  observed  at  a  tropopanse  height  of  ii  =  6.5  km.  For  the  nnmerical  simnlations,  a 
higher  tropopanse  height  will  take  longer  to  reach  steady-state  than  a  lower  height,  so  the 
tropopanse  height  of  ii  =  6.5  km  was  selected  for  onr  nnmerical  investigation  of  tropopanse 
tnning. 

In  the  nnmeric  code,  the  Brnnt-Vaisala  freqnency  prohle  is  slightly  smoothed.  The  fnnctional 
form  is  a  hyperbolic  tangent  specihed  as 

N  =  ^{3  +  ttN  +  bNtdJih[cN  {z  -  cIn)]}  (113) 

where  the  fonr  parameters  a^,  b^,  cjq  and  determine  the  location  of  the  tropopanse  and 
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Figure  22:  A  comparison  of  numerically  calculated  nonlinear  streamlines  with  the  linear,  analytic  solution 
for  two  layers  of  uniformly  stratified  topography.  The  tropopause  is  located  at  a  non-dimensional  height 
oi  H  =  6.5.  In  the  numeric  code,  this  is  done  using  the  functional  form  of  (113).  The  amplitude  of  the 
topography  is  e  =  =  0.4  and  the  Rossby  number  is  Ro  —  10. 


the  sharpness  of  the  transition  of  Brnnt-Vaisala  freqnency  in  going  from  the  troposphere  to 
the  stratosphere. 

Figure  22  shows  the  streamlines  for  the  numeric  solution  of  stratihed  flow  over  topography 
in  the  presence  of  weak  rotation,  and  provides  a  comparison  with  the  linear,  analytic,  hydro¬ 
static  two-layer  streamfunction  solution  of  (87).  The  linear,  analytic  and  nonlinear,  numeric 
solutions  differ  considerably.  In  addition  to  numeric  streamlines  being  shifted  downstream  of 
the  linear,  analytic  streamlines  (which  is  partially  due  to  non-hydrostatic  effects),  the  non¬ 
linear  slopes  are  also  steeper.  Because  the  numerical  results  are  limited  by  computational 
constraints,  it  is  impossible  to  see  from  these  results  whether  the  inertial-gravity  waves  prop¬ 
agate  far  downstream  as  they  do  in  the  analytical  solution. 

Figure  23  shows  the  spanwise  velocity  calculated  for  two  layers  of  uniformly  stratihed  how 
with  a  tropopause  height  oi  H  =  6.5  km.  The  top  hgure  shows  the  linear,  hydrostatic,  ana¬ 
lytic  response  and  the  lower  hgure  shows  the  fully  numeric,  nonlinear  solution.  Qualitatively 
the  results  show  that  the  nonlinear  response  in  the  troposphere  is  much  stronger  than  the 
corresponding  linear  response.  This  is  consistent  with  the  large  nonlinear  forcing  coefficients 
in  the  troposphere  of  the  analytic  solution. 

Figure  24  provides  a  more  quantitative  comparison  of  the  spanwise  velocity  for  the  linear  an 
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Figure  23:  Spanwise  velocity  profile  obtained  from  (a)  the  linear,  analytic  solution  (b)  the  nonlinear,  numeric 
solution  for  parameters  given  in  Figure  22. 
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Figure  24:  Spanwise  velocity  as  a  function  of  altitude  for  several  locations  downstream  of  of  the  topography 
for  parameters  given  in  Figure  22.  Dimensions  of  downstream  locations  are  in  kilometers,  i.e  40  km,  50  km, 
etc. 

nonlinear  cases.  In  the  figure,  the  vertical  profile  of  spanwise  velocity  is  plotted  at  several 
stations  downstream  of  the  topography.  Above  the  tropopause,  there  is  very  little  difference 
between  the  linear  and  nonlinear  solutions.  However,  below  the  troposphere,  the  nonlinear 
spanwise  velocity  is  much  larger  (in  magnitude)  than  the  linear  response,  which  is  consistent 
with  the  qualitative  observations  made  in  Figure  23. 

In  order  to  make  an  even  more  detailed  comparison  between  the  linear  and  nonlinear  so¬ 
lutions,  the  vertical  spectrum  of  the  spanwise  velocity  for  both  solutions  was  computed  at 
several  locations  downstream  of  the  topography.  The  results  are  shown  in  Figure  25.  In 
the  lower  layer,  the  propagating  component  of  the  linear  solution  should  have  a  vertical 
wavenumber  of  unity.  For  this  case,  the  vertically  propagating  component  of  the  linear  solu¬ 
tion  in  the  stratosphere  should  have  a  vertical  wavenumber  of  two  because  the  Brunt-Vaisala 
frequency  in  the  stratosphere  is  twice  that  in  the  troposphere.  The  linear  spectra  have  strong 
peaks  at  the  appropriate  wavenumbers  and  are  essentially  zero  for  higher  harmonics,  which  is 
consistent  with  the  formulation  of  the  linear  solution  given  in  (87).  The  forcing  of  the  outer 
solution  by  nonlinearity  in  the  inner  solution  gives  rise  to  higher  harmonics  in  the  vertical, 
as  formulated  in  (108).  These  higher  harmonics  are  expressed  using  the  Fourier  coefficients 
(0)  and  (0).  The  vertical  spectra  of  spanwise  velocity  downstream  of  the  topography 
in  the  troposphere  show  these  higher  harmonics.  However,  the  spectra  in  the  stratosphere 
do  not.  Again  this  is  consistent  with  Figure  24  and  indicates  that  the  nonlinearity  in  the 
inner  solution  gives  rise  to  higher  harmonics  in  the  troposphere  that  propagate  far  down- 
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stream,  but  the  response  in  the  stratosphere  is  largely  unaffected  by  the  nonlinearity  for  this 
particnlar  tropopause  height.  A  comprehensive  stndy  of  the  effects  of  tropopause  altitnde 
has  not  been  done.  It  is  possible  that  snch  a  study  would  reveal  tropopause  heights  giving 
rise  to  more  pronounced  nonlinearity  in  the  stratosphere  and  less  in  the  troposphere. 
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Figure  25:  Vertical  spectrum  of  spanwise  velocity  at  several  horizontal  positions  downstream  of  the  to¬ 
pography  for  parameters  given  in  Figure  22.  Horizontal  positions  downstream  of  topography  are  given  in 
kilometers,  i.e  40  km,  50  km,  etc.  (a)  Spectrum  in  the  stratosphere,  (b)  Spectrum  in  the  troposphere. 
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3  Effect  of  Unsteady  Wind 
3.1  Introduction 

The  classic  analyses  of  Queney  [17]  and  Long  [14]  assume  that  flow  of  stratihed  fluid  over  to¬ 
pography  is  steady  in  time.  Many  of  the  subsequent  studies  have  also  used  this  assumption, 
and  there  are  only  a  handful  of  investigations  that  account  for  unsteady  wind  in  strati¬ 
hed  how  over  topography  in  the  atmosphere.  Among  these  are  the  theoretical  works  of 
Bell  [4] ,  Bannon  and  Zehnder  [1]  and  Hines  [9]  and  the  numerical  investigation  of  Lott  and 
Teitelbaum[15].  More  recently  the  ehect  of  unsteady  wind  over  topography  has  been  exam¬ 
ined  by  Skopovi[19].  Based  on  fully  numerical  simulations,  Skopovi  pointed  out  that  slowly 
accelerated  transient  wind  could  produce  streamlines  that  were  signihcantly  steeper  than 
those  for  an  impulsively  started  wind.  In  addition,  his  work  showed  that  a  monochromatic 
temporal  huctuation  in  wind  speed  could  produce  waves  with  larger  amplitudes  than  those 
for  a  steady  wind  corresponding  to  the  mean  wind  speed.  Skopovi  [19]  used  a  fully  numerical 
code  to  examine  the  ehect  of  an  oscillatory  wind  component, 

U  =  Uq  (1  -|-  A  cos  ut) , 

on  nonlinear  mountain  waves.  He  noted  that,  for  low  background  how  frequencies,  uj/Nq  = 
0.1,  say,  (Ao  is  the  buoyancy  frequency),  the  nonlinear  unsteady  response  may  be  amplihed 
and  can  be  dramatically  diherent  from  the  steady-state  response  corresponding  to  the  mean 
wind  Uq.  This  is  illustrated  in  Fig.  26  below,  which  shows  how  streamlines  at  t  ~  0.5  hr  and 
t  ~  2.2  hr  after  the  wind  has  been  turned  on.  Note  that,  while  the  wind  speed  is  equal  to 
Uq  =  10  m/sec  at  both  these  times,  the  response  at  t  ~  2.2  hr  is  signihcantly  stronger  and 
quite  diherent  from  the  steady-state  response. 

These  results  call  for  a  systematic  study  of  the  question:  why  does  slowly  varying  wind 
matter?  Qualitatively,  the  answer  can  be  traced  to  the  fact  that  the  response  in  the  hydro¬ 
static  limit  (/i  =  Uq/  {NL)  1)  evolves  on  a  relatively  long  timescale  t  ~  As  a 

result,  low-frequency  wind  oscillations  can  have  an  O  (1)  ehect,  particularly  in  the  nonlinear 
response.  Moreover,  as  the  group  velocity  vanishes  in  the  frame  of  the  mountain,  the  hydro¬ 
static  response  is  resonant,  so  small-amplitude  wind  oscillations,  A  =  can  have  an 

O  {1)  ehect  as  well.  In  this  investigation  we  present  an  analytic  theory  for  stratihed  how  over 
topography  in  the  presence  of  an  unsteady  background  how  along  with  some  fully  numeric 
computational  results. 


3.2  Governing  Equations 

The  asymptotic  theory  of  Kantzios  and  Akylas  [10],  hereafter  referred  to  as  KA,  generalizes 
Long’s  solution  to  account  for  slow  temporal  variations  and  weak  spatial  variations  in  up¬ 
stream  wind  speed  and  Brunt-Vaisala  frequency.  The  theory  of  KA  is  analogous  to  that  of 
Grimshaw  and  Yi  [8],  who  examined  the  case  of  a  hnite  depth  how  over  topography.  KA 
gives  the  leading-order  streamfunction  disturbance  as 

'4’  =  y  +  -i-  c.c.^  (114) 
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Figure  26:  Previous  numerical  results  of  Skopovi  [19]  depicting  response  streamlines  due  to  slowly  fluctuating 
wind  of  the  form  U  =  [/q  (1  +  Acoswt),  where  Uq  =  10  m/s,  A  =  0.5  and  uj/Nq  =  0.1,  over  topography  of 
peak  height  h  =  0.6  km.  (a)  Fluctuation  of  free-stream  flow  in  time.  Results  in  Fig.  2.b  are  both  at  times 
when  the  free-stream  flow  is  at  10  m/s  but  are  separated  by  one  period  of  fluctuating  free-stream  velocity, 
(b)  Comparison  of  transient  response  at  two  different  times,  which  both  correspond  to  a  free-stream  velocity 
of  10  m/s.  Despite  the  free-stream  velocity  being  the  same  in  both  cases,  the  streamlines  differ  significantly, 
indicating  the  importance  of  the  fluctuating  component  of  the  wind,  (c)  A  comparison  of  the  transient 
response  for  a  fluctuating  wind  to  the  steady-state  response  for  a  steady  wind  given  by  Long’s  solution. 
The  results  show  that  wave  amplitudes  forced  by  the  transient  wind  can  be  significantly  larger  than  those 
obtained  in  the  steady-state  solution  for  a  steady  wind.  This  unsteady  wind  may  lead  to  wave  breaking  at 
mountain  heights  smaller  than  those  predicted  by  Long’s  steady-state  solution. 
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where  A  is  the  complex  amplitude  function 


A  =  a(x,  Y,  T)  +  ib  {x,  Y,  T) 


(115) 


and  time  and  altitude  have  been  scaled  with  the  square  of  the  longwave  parameter,  p  = 
<  1  as 

T  = 


Y  =  fi^y 


(116) 


where  y  and  t  are  the  nondimensional  altitude  and  time  dehned  previously.  Then  the  complex 
amplitude  is  governed  by  a  system  of  Volterra  integral  equations.  Upon  taking  one  derivative 
with  respect  to  x  of  the  hrst-kind  formulation  given  in  KA,  the  evolution  equations  become 
Volterra  equations  of  the  second  kind  given  by 


KiittT  +  A  f  {^iixO-T  +  Ki2xb''j') 

J  —OO 


1  2 

Ojxxx  T  bxY  0 


/X 

dx  {K21x(It  Y  K22xb'j’) 

-OO 


—  a-xY  =  0 


(117) 


where  the  kernels,  Kij  are  in  general  nonlinear  functions  of  the  complex  amplitude,  A  =  a+ib, 
and  are  given  in  Appendix  F.  The  evolution  equations  are  subject  to  the  kinematic  boundary 
condition  at  the  surface  of  the  topography  requiring  that  the  curve  of  topography  be  a 
streamfunction.  Doing  so  gives  the  boundary  condition 

acose/ —  6sine/ =  — -e/  (118) 

where  y  =  ef  {x,T)  is  the  curve  dehning  the  topography. 


Following  the  work  of  Bell  [4] ,  we  transform  the  governing  equation  from  one  in  which  the 
mountain  is  hxed  and  the  wind  unsteady  to  one  in  which  the  wind  is  steady  and  the  mountain 
oscillates.  Assuming  a  spatially  uniform,  undisturbed  wind  of  the  form 

V  =  Uo(l  +  /i^AcosDT)  (119) 

where  Uo  is  the  temporally  averaged  wind  speed,  A  is  an  amplitude  of  order  unity  and  D  the 
frequency  of  order  unity.  Because  p  -C  1,  the  wind  is  nearly  constant  in  time  with  a  very 
small  fluctuating  component.  The  frequency  D  multiplies  the  stretched  time,  T,  so  despite 
D  =  (9  (1),  the  frequency  in  unstretched  coordinates  is  very  small  (p^D).  The  combination 
of  a  small  amplitude  wind  oscillation  at  low  frequency  can  give  a  leading-order  disturbance. 
In  the  transformed  coordinates  where  the  wind  is  steady,  this  corresponds  to  the  mountain 
having  an  amplitude  of  oscillation  of  order  unity. 


We  dehne  a  new  horizontal  coordinate,  for  which  the  wind  is  steady  and  the  topography 
oscillates.  For  an  unsteady  wind  of  the  form  (119),  the  transformation  is 


X  + 


A 

Q 


sinDT 


(120) 
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where  the  amplitude  of  the  mountain  oscillations  in  the  steady  wind  coordinate  system  is 

—  =  (9  (1),  which  leads  to  a  large-amplitude  unsteady  response  despite  the  small-amplitude 

oscillations  of  wind  in  the  original  coordinate  system.  In  this  new  coordinate  system,  the 
curve  dehning  the  topography  becomes 


y  =  ef  {x,T)  =  ef  (0  =  e 


(Gaussian) 


1 

r+e’ 


(Witch  of  Agnesi) 


(121) 


3.3  Linear  Limit 

In  the  limit  as  the  amplitude  of  the  topography  goes  to  zero,  e  — 0,  the  matrix  of  kernels  in 
(117)  reduces  to  the  diagonal  matrix  as 

All,  A22  — ^  1,  Ki2,  K21  — 0  as  e  — 0  (122) 


In  this  limit,  the  evolution  equations  are 
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(It  ~  ~  ^ 


br  —  —  a^y  —  0 


}  1  e 


(123) 


subject  to  the  linearized  boundary  condition 


a  =  --e/(0,  l"  =  0,e^0 


(124) 


If  we  further  make  the  assumption  that  the  response  is  hydrostatic,  then  the  dispersive  terms, 
and  b^^^  become  zero  and  complex  amplitude.  A,  can  be  written  as  an  inverse  Fourier 
transform  as 


A  =  —e 
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f{k)  e^'^^exp 
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dk,  fi,e  ^  0 


(125) 


Skopovi  [19]  has  also  derived  an  alternate  linear  solution  for  the  non-hydrostatic,  steady-in- 
the-mean  flow  over  topography.  The  solution  incorporates  dispersive  effects  and  neglects  the 
initial  wave  propagation  in  the  transition  from  an  undisturbed  flow  to  the  oscillatory  flow 
over  topography.  Rather,  the  solution  gives  the  oscillatory  flow  over  topography  once  it  has 
reached  a  periodic  state.  The  solution  is 
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where  Uq  is  the  frequency  in  the  unstretched  time  domain,  t,  and  the  vertical  wavenumber, 
m,  is  chosen  to  satisfy  boundedness  and  the  radiation  condition,  and  is  given  as 
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3.4  Nonlinear  Asymptotic  Solution 

For  the  general  case  of  a  hnite-amplitude  topography,  the  matrix  of  kernels  does  not  reduce 
to  the  diagonal  matrix.  For  weakly  nonlinear  solutions,  the  kernels  can  be  evaluated  analyt¬ 
ically  through  asymptotic  expansion  of  (114).  However,  at  the  limit  of  wave  breaking  when 
the  streamlines  become  vertical,  the  kernels  become  singular.  As  this  limit  is  approached, 
the  kernels  must  be  evaluated  numerically.  In  addition  to  the  presence  of  the  off-diagonal 
kernels  and  the  x— derivatives  of  the  kernels  in  the  integrals,  the  fully  nonlinear  problem 
differs  from  the  linear  problem  in  that  the  boundary  condition  must  be  applied  along  the 
contour  of  the  topography  as  in  (118),  rather  than  at  y  =  0,  as  in  (124).  Several  numerical 
methods  have  been  investigated  for  the  integration  of  these  nonlinear  evolution  equations. 
Most  notable  among  them  was  the  modihed  Lax-Wendroff  method  employed  by  Prasad  and 
Akylas  [16].  However,  in  each  of  the  methods  tried,  grid-scale  oscillations  appeared  for  to¬ 
pography  of  amplitude  greater  than  approximately  e  =  0.50.  These  grid-scale  oscillations 
generally  appeared  after  approximately  one  quarter  of  the  period  of  oscillation  of  the  wind, 
which  was  on  the  order  of  one  hundred  computational  timesteps.  These  grid-scale  oscillations 
were  not  seen  for  smaller-amplitude  topography.  The  method  of  hltering  used  by  Prasad 
and  Akylas  [16]  was  implemented  in  an  attempt  to  eliminate  these  grid-scale  fluctuations. 
However,  mild  hltering  was  unable  to  adequately  eliminate  these  erroneous  oscillations,  and 
work  on  their  elimination  for  larger-amplitnde  topography  is  ongoing. 


Figure  27  shows  a  streamlines  for  unsteady  how  over  topography  for  a  nondimensional  monn- 
tain  amplitnde  of  e  =  0.5  and  amplitnde  of  wind  oscillation  that  is  18  percent  of  the  mean 
wind  speed  of  10  m/s.  A  comparison  of  hve  diherent  solntions  is  presented.  The  solntion 
labeled  “Linear  Long’s”  is  Long’s  solution  for  steady  wind  in  the  linear  limit.  The  solntion 
labeled  “Long’s”  is  Long’s  solution  for  steady  wind  over  hnite-amplitude  topography.  The 
solution  labeled  “CFD”  is  the  fully  numeric  solution  obtained  by  integrating  the  governing 
eqnations  for  the  primitive  variables  {u,  v  p  and  p)  in  time.  The  solution  labeled  “Asymp¬ 
totics”  is  the  integration  of  the  governing  eqnations  given  in  (117)  snbject  to  the  bonndary 
condition  (118).  Finally,  the  solntion  labeled  “Linear  Asymptotics”  is  the  solution  to  the 
linearized  governing  eqnations  in  (123)  snbject  to  the  linearized  bonndary  condition  (124). 

Above  the  topography,  the  agreement  between  the  asymptotic  theory  of  KA  and  the  com- 
pntational  solution  is  reasonably  good.  Downstream  and  aloft,  the  fully  numeric  solution 
exhibits  waves  of  wavelength  approximately  12  km.  For  this  relatively  small-amplitnde  to¬ 
pography  (e  =  0.50),  the  nnsteady  streamlines  do  not  appear  to  be  signihcantly  steeper  than 
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Figure  27:  Streamlines  for  unsteady  flow  over  topography.  ^  =  g-  e  =  0.5.  A  =  Q  —  6.7858.  Computational 
time  t computational  —  250. 
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those  for  the  steady  flow  of  Long’s  solution.  However,  they  are  considerably  steeper  than 
the  corresponding  linear  streamlines,  suggesting  that  as  the  amplitude  of  the  topography  in¬ 
creases,  the  nonlinear  effect  of  unsteadiness  increases.  Extrapolating  this,  it  is  plausible  that 
slightly  unsteady  wind  can  produce  signihcantly  steeper  streamlines,  as  shown  by  Skopovi 
[19]. 

Despite  not  having  an  asymptotic  solution  for  larger  amplitudes,  numerical  simulations  have 
been  done  for  unsteady  wind  over  mountains  with  amplitudes  larger  than  e  =  0.50.  Figure 
28  shows  streamlines  for  the  fully  numeric  simulation  of  stratihed  flow  over  topography  with 
unsteady  upstream  wind.  The  parameters  are  the  same  as  those  in  the  previous  case  with  the 
exception  that  the  mountain  height  is  now  ho  =  0.73  km,  corresponding  to  a  nondimensional 
amplitude  of  e  =  0.73  for  characteristic  wind  speed  Uq  =  10  m/s  and  Brunt-Vaisala  frequency 
Nq  =  0.01  sec“^.  The  results  show  that  the  nonlinear  response  with  weak  fluctuations  in 
upstream  wind  speed  is  signihcantly  stronger  than  that  predicted  by  Long’s  solution  for 
steady  wind.  In  particular,  the  streamline  at  approximately  a;  =  10  km,  y  =  A  km  is 
quite  steep  compared  to  Long’s  solution.  This  suggests  that  the  critical  amplitude  for  wave 
breaking  for  the  unsteady  wind  is  smaller  than  the  critical  amplitude  for  wave  breaking  for 
the  steady  Long’s  solution. 
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Figure  28:  Streamlines  for  unsteady  flow  over  topography.  fJ-  =  ^-  e  =  0.73.  A  —  Q  =  6.7858.  Computational 
time  t computational  —  450. 
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4  Conclusions 


•  In  stratified  flow  over  topography,  rotation  behaves  as  a  singular  perturbation  when 
Ro  3>  1.  For  mountains  of  width  less  than  approximately  L  =  50  km,  Coriolis  forces  do 
not  have  sufficient  time  to  act  on  the  flow  as  it  passes  over  the  mountain.  However,  far 
downstream  of  the  mountain,  rotation  becomes  as  important  as  stratihcation  because 
the  Coriolis  force  has  sufficient  time  to  act  upon  the  flow.  The  results  of  §1  show  the 
effect  of  rotation  far  downstream  of  topography  for  large  Rossby  number. 

•  A  matched  asymptotic  expansion  has  been  developed  for  nonlinear  stratihed  flow  over 
topography  in  the  limit  Ro  3>  1.  The  inner  region  of  the  domain  is  near  the  mountain 
horizontally  and  extending  upward  indehnitely.  In  this  region,  rotation  is  a  weak  per¬ 
turbation  and  the  leading-order  solution  is  that  of  Long.  The  outer  region  of  the  domain 
is  far  from  the  topography:  a;  3>  Here  the  solution  is  linear  but  the  rotation,  in 
addition  to  stratihcation,  is  of  leading-order  importance.  The  matching  process  shows 
that  nonlinear  interaction  in  the  inner  solution  give  rise  to  higher-harmonic  oscillations 
that  propagate  far  downstream  of  the  topography. 

•  Nonlinear  interactions  above  the  mountain  give  rise  to  steady,  higher-harmonic  gravity- 
inertial  waves  far  downstream.  The  effect  of  nonlinearity  for  a  single  layer  is  most  clearly 
seen  in  comparing  the  linear  spanwise  velocity  shown  in  Figure  8  with  the  corresponding 
nonlinear  spanwise  velocity  depicted  in  Figure  9.  The  structure  of  the  spanwise  velocity 
for  the  nonlinear  case  shows  closed  contours  indicative  of  higher  harmonics  in  the  ver¬ 
tical  direction.  These  higher  harmonics  are  the  signature  of  the  nonlinear  interactions 
above  the  topography.  For  the  single  layer,  these  nonlinear  interactions  are  relatively 
weak  because  Long’s  solution  (with  the  nonlinear  boundary  condition)  does  not  differ 
appreciably  from  the  linear  limit  of  Long’s  solution.  Thus  nonlinear  interactions  are 
relatively  weak  and  the  higher  harmonics  far  downstream  are  relatively  weak  for  the 
single-layer  model  of  the  atmosphere.  The  effect  becomes  stronger  when  the  tropopause 
is  taken  into  account. 

•  The  asymptotic  theory  has  been  qualitatively  validated  via  comparison  with  fully  nu¬ 
meric  simulations  of  stratihed  how  on  an  /—plane.  In  particular,  the  simulations  for 
two  layers  of  uniformly  stratihed  how  show  the  presence  of  higher  harmonics.  This 
can  be  seen  in  the  vertical  spectrum  of  spanwise  velocity  in  Figure  25.  However,  com¬ 
parison  of  the  analytic  theory  with  these  numeric  simulations  is  imperfect.  In  order 
to  compare  with  the  asymptotic  theory,  the  numeric  simulations  should  give  steady- 
state  solutions  far  downstream  of  the  topography.  Because  of  limits  in  computational 
resources,  a  steady-state  solution  is  only  obtained  for  relatively  short  distances  down¬ 
stream  {x  <  900  km),  preventing  a  good  comparison  with  the  analytical  results.  Never¬ 
theless,  the  higher  harmonics  in  the  spanwise  velocity  shown  in  Figure  25  suggest  that 
nonlinear  interactions  above  the  topography  are  signihcant  and  do  give  rise  to  higher 
harmonics  downstream  of  the  topography. 

Because  the  analytical  results  can  be  easily  evaluated  arbitrarily  far  downstream,  the 
analytical  theory  provides  a  predictive  model  that  is  valuable  when  numerical  simula¬ 
tions  may  be  impractical  due  to  computational  limitations. 
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•  The  combined  effect  of  the  Earth’s  rotation  and  the  tropopause  is  to  enhance  and 
extend  the  wave  activity  far  downstream  from  the  topography,  as  shown  in  §2.  The 
effect  of  tropopause  tuning  for  two  different  tropopause  heights  was  examined  in  §2. 
For  a  tropopause  height  oi  H  =  12.5  km  (based  on  a  wind  speed  Uq  =  10  m/s  and 
Brunt-Vaisala  frequencies  of  Ni  =  0.01  sec“^  and  N2  =  0.02  sec“^  in  the  troposphere 
and  stratosphere,  respectively)  the  response  above  the  mountain  had  signihcantly  larger 
streamline  slopes  and  nonlinear  interactions  than  the  case  with  a  tropopause  height  of 
H  =  11.0  km.  These  nonlinear  interactions  above  the  topography  produced  higher 
harmonics  in  the  wavefield  far  downstream  of  the  topography,  as  shown  in  Figures  18 
and  19.  In  general,  the  tropopause  height  at  which  the  response  is  tuned  depends  on 
the  Brunt-Vaisala  frequencies  in  the  troposphere  and  stratosphere  as  well  as  on  the 
velocity  upstream  of  the  topography. 

•  The  response  due  to  unsteady  wind  can  be  signihcantly  stronger  than  the  response 
due  to  steady  wind.  The  analytic  theory  shows  that,  even  for  small  huctuations  in 
wind  speed,  waves  of  signihcantly  larger  amplitude  can  be  produced  than  those  for  a 
steady  wind  corresponding  to  the  mean  wind  speed.  The  amplitude  of  the  response  is 
based  on  the  ratio  of  the  amplitude  of  huctuation  in  wind  speed.  A,  to  the  frequency 
of  the  huctuation,  11,  as  given  in  Eqn.  (119).  Thus  for  low-frequency,  low-amplitude 
huctuations,  the  unsteady  response  can  be  of  order  unity.  This  unsteady  wind  can  cause 
wave  breaking  at  lower  topography  steepness.  The  previous  numeric  results  of  Skopovi 
26  as  well  as  the  numeric  results  of  the  current  study  shown  in  Figures  27  and  28  conhrm 
that  an  unsteady  wind  can  produce  steeper  streamlines  than  the  corresponding  steady 
wind. 

Future  work  includes  the  following: 

•  The  preceding  analysis  of  tropopause  tuning  only  considered  the  comparison  of  two 
tropopause  heights.  The  tuned  case  yielded  a  signihcantly  amplihed  response  in  the 
troposphere,  but  the  response  in  the  stratosphere  was  not  signihcantly  enhanced  by  the 
nonlinearities.  It  is  believed  that  other  tropopause  heights  will  yield  larger  nonlinear¬ 
ities,  and  a  correspondingly  larger-amplitude  response  in  the  stratosphere,  and  a  less 
signihcant  amplihcation  in  the  troposphere.  A  search  for  tropopause  heights  producing 
amplihed  waves  in  the  stratosphere,  as  opposed  to  the  case  shown  in  Figure  19  for  which 
amplihcation  occurs  in  the  troposphere,  will  be  done. 

•  Thus  far,  asymptotic  results  for  unsteady  wind  for  topographies  of  amplitude  e  >  0.5 
have  not  been  obtained.  The  asymptotic  model  discussed  in  §3  will  be  used  to  study 
the  ehect  of  unsteady  wind  over  steeper  topographies.  These  results  will  be  compared 
with  the  fully  numeric  simulations. 

•  The  asymptotic  model  discussed  in  §3  is  valid  for  slowly  oscillating  wind  with  small- 
amplitude  oscillations.  The  current  model  will  be  extended  to  relax  this  restriction  of 
small-amplitude  wind  oscillations,  and  a  study  of  large-amplitude  wind  oscillations  over 
hnite-amplitude  topography  will  be  made. 
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III.  APPENDICES 


A  Inner  Limit  of  Outer  Solution 

The  outer  solution  is 
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Defining  the  four  integrals  in  the  solution  as 
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Taking  the  limit  as  5;  — >  0  gives 
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Using  the  theory  of  generalized  functions  [13],  we  rewrite  these  integrals  as 
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In  hnding  the  limit  of  1^  as  x  ^  0,  the  same  transformation  is  made  as  for  12-  Then 
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where  Kq  is  the  modihed  Bessel  function  of  the  second  kind.  The  asymptotic  behavior  of 
the  modified  Bessel  function  of  the  second  kind  is 
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Taking  the  limit  of  the  second  integral  in  /4  as  5;  — 0  gives 
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The  sgna;  function  results  in  a  jump  in  the  spanwise  velocity  w  across  x  =  0,  giving 
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B  Matching:  Jump  in  Spanwise  Velocity  due  to  Inner  Solution 


The  leading  order  spanwise  velocity  in  the  inner  solution  is 
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We  would  like  to  determine  the  jump  in  spanwise  velocity  in  going  from  x— >— cx)tox— >cx). 
This  jump  will  be  necessary  in  order  to  match  the  inner  solution  with  the  outer  solution. 
We  start  by  analyzing  the  behavior  of  the  integral  containing  h{x).  From  the  radiation 
condition,  b  (x)  =  {a  (x)}  where  is  the  Hilbert  transform.  Then 
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The  hrst  integrand  is  odd  in  x  and  the  second  term  even  in  x.  Looking  at  the  integral  from 
—oo  to  cx),  the  even  integrand  vanishes  and  we  are  left  with 
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because  each  of  the  conjugates  is  purely  imaginary.  Thus  the  term  containing  b  (x)  does  not 
contribute  to  a  finite  jump  across  the  topography  in  going  from  x— cx)tox— s^cx).  Then 
we  are  left  with  a  hnite  jump  in  the  spanwise  velocity  given  by 
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C  Details  of  Matching  -  Single  Layer 


In  order  to  match  the  outer  solution  to  the  inner  solution,  we  first  scale  the  leading  order 
outer  solution  by  which  is  unity  in  the  linear  limit  where  a  (x)  =  —  \eh{x).  Then  the 

outer  solution  becomes 
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with  asymptotic  behavior 
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The  In/?  term  in  the  expression  for  w  can  be  made  to  agree  with  the  inner  solution  by 
dehning  the  constant  of  integration  in  the  inner  solution  as 
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At  this  point,  the  outer  solution  does  not  yet  satisfy  the  required  jump  in  spanwise  velocity 
across  the  inner  solution.  This  jump  is  induced  by  the  higher  harmonics,  [x)  and  (x). 
Thus  we  must  add  another  component  to  the  outer  solution  to  achieve  the  necessary  jump 
in  spanwise  velocity.  We  write 
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Since  the  outer  problem  is  linear,  we  can  superimpose  solutions  and  separate  and  wn 
into  vertical  Fourier  components  as 
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such  that  each  harmonic  satisfies  the  appropriate  matching  conditions  with  the  inner  solution, 
namely 
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A  new  variable,  can  be  defined  such  that 
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=  — - —  [fn  (0)  sin  ny  +  (0)  cos  ny]  +  F  {k)  n  >  1 


k^-kl 


where 


kn.  — 


n 


n^  —  1 


>  1 
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and  the  vertical  wavenumber  is 


m  {k)  = 


k 


:  \k\ 

'y/l-k^' 


(1^1  <1) 


which  is  determined  by  the  radiation  condition.  The  boundary  condition  of  zero  streamline 
disturbance  at  y  =  0  requires 


Therefore,  the  solution  is 


^  Inf  r  n  ^ 

[f-n  (0)  sin  ny  +  Sn  (0)  [cos  ny  -  }  ,  n  >  1 


(137) 


In  inverting  the  Fourier  transform,  care  must  be  exercised  because  has  poles  on  the  real 
axis  at  k  =  tzkn-  The  contours  must  be  deformed  in  a  manner  consistent  with  the  radiation 
condition  in  the  horizontal  direction.  Since  steady  waves  are  observed  downstream  of  the 
topography  only,  we  require 


^9  X 


\k\  1 

m? 


Thus  the  integral  is  deformed  beneath  the  poles  giving 

=  — 27r  cos  knX  [fn  (0)  sin  ny  +  Sn  (0)  cos  ny]  H  [x) 


Dehning 


+rsn  (0) 


oo  u^ikx 


k‘^-kl 


exp  i 


dk>,  n  >  1 


4”^  =  s 


’1  U„ikx 


dk\  ,  n  >  1 


jin)  —  fx, 
Jg  — 


we  rewrite  the  solution  as 


oo  i.Akx 


U2  _  U2 
1  ^ 


exp  i 


dk  \  ^  n>  1 


=  — 27rcos  knX  [r„  (0)  sinny  +  (0)  cosny]  H  (x)  +  2s„  (0)  ([4"’^  +  7^^)  ,  n  >  1 


and  examine  the  limit  as  5;  — 0.  For  the  hrst  integral 


lim  D"'  = 


k 

,  k^-kl 


exp  — 


dk}  =  t) 
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because  the  integrand  converges  to  a  real  number  as  x  — 0.  Next  we  make  the  same 
transformation  for  as  was  made  in  §1.4  where  a  =  \/k?‘  —  1  \x\.  Then 


oo  ^i\/ cr2+a;2  gg^  x 


0  -  (^n  -  1) 


i\^a^  + . 


As  5;  — 0,  the  integral  becomes 


ly  ^^lexp  (iy) 


roc  (je*®"  ^ 

0  a2-(A;2-l)p‘ 


da  ?  ,  5;  — >  0 


Now  the  integral  behaves  differently  depending  on  whether  x  is  greater  or  less  than  zero.  For 
x  >  0  we  close  the  contour  of  integration  to  the  top  and  there  is  a  pole  inside  the  contour, 
giving 

~  ^  jexp  (iy)  m+  [  .  ^ I  ’  ^  ^ 

I  L  Jo  p^  +  {,kl-l)x^  \] 

For  negative  x,  the  contour  can  be  closed  to  the  bottom,  thereby  excluding  the  singularity 
at  a  =  J y  —  lx  and  giving 


-dpi  [  ,  X  ^  0^ 


ly  ~  A  exp  {iy) 


.g  p2  +  ^kl-l)x^^y^ 

For  both  5;  >  0  and  5;  <  0  it  is  then  necessary  to  evaluate  the  integral 


L  p^  +  {y  - 1) x^"^^  2r 


A\/k^—lxi 


Eii  [i\  kl^-  lx]  +  e 


—ix/kl  —  lxi 


Eii  —  lx 


=  -7-ln  (^.Jy-l\x\^  +  ^7r^Jy-lx  +  0{(kl  -  l)  5^}  , 


Therefore 


ly  ~  A  {exp  {iy)  [mH  {x)  —  In  |5;|]} 


~  nH  {x)  cosy  —  In  \x\  siny 
Finally,  the  behavior  of  as  5;  — 0  is 

~  — 27r  [fn  (0)  sin  ny  +  Sn  (0)  cos  ny]  H  {x) 


+2sn  (0)  [Trif  {x)  cosy  —  In  \x\  siny] ,  |5;|  — 0,  n  >  1 
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Taking  derivatives  with  respect  to  x  and  y  gives  and  w,  respectively. 

^(n) _ 

X 

~  —27171  [fn  (0)  cos  ny  —  Sn  (0)  sin  ny]  H  {x) 

—2sn  (0)  [tiH  (x)  siny  +  In  \x\  cosy] 

=  27771  [—Tn  (0)  COS uy  +  Sn  (0)  sin ny]  —  27rs„  (0)  siny,  n  >  1 
This  result  contains  the  necessary  jump  in  spanwise  velocity  across  the  topography  of 

wIq-  =  27171  [—fn  (0)  cosny  +  Sn  (0)  sinny] 

but  it  also  contains  an  additional,  unwanted  term  in  the  jump.  Furthermore,  w  and  ip  do 
not  decay  sufficiently  fast  as  5;  — 0.  Thus  in  order  to  eliminate  these  additional,  unwanted 
terms,  we  add  a  linear  multiple  of  the  original  linear  solution,  which  was  previously  given  as 


The  solution  for  the  first  harmonic  outer  solution  must  satisfy 

(</■"  +  </■«),,  +  =  0 

=  0,  (y  =  0,5;  7^0)  (139) 
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subject  to  the  matching  conditions  as  5;  — 0  given  by 


=  o  (In  |a;|)  J 


w 


(1) 


0+ 


We  then  write  the  solution  as 


_  =  271  [ri  (0)  cos  y  —  §1  (0)  sin  y] 


^(1)  =  +  ^(1) 


(1)  _  ,7,(1)  ^  ,7,(1) 


where 


b'’  =  S 


’  =  wy  + 


exp  ( \  e^'^^dk  +  f  cos 


ky 


wW  = 


'•I  ^ikx 


=  ^\  sin 
-ky  \ 


ky  \  ^iki 


—  1 

dk  + 


oo  ^ikx 


:  sin 


ky 


~ (1) 

wl  =  ^ 


^oc  ^ikx 


COS 


1  y/y-1 

^  ky 


dk 


such  that  both  and  independently  satisfy  (139).  It  is  now  necessary 

to  determine  the  asymptotic  behavior  of  this  solution  as  5;  — 0.  We  dehne 


I7  =  ?R: 


h  =  -^ 


cos 


ky 


'•OO  ^ikx 


:  Sin 


ky 


such  that 


V'b  =  /,  +  /, 


dk 


ini')  =  h  +  h 

which  utilizes  Ji  and  Is  dehned  in  Appendix  A.  Then  making  the  transformation  a  = 
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—  1 1 5;  I  as  before  gives 

r  ^  roo 


17-3? 

Jg  ~  — ^  |sin|/ 


X 


oo  ^icF  sgn  X 


cosy 

\x\ 


^0-2  + 


da  >  =  —  sin  y  sgn  5; 


cos  ada  rsj  Q 


sin  a  TT 

d<7  ——  sgn  a;  sin  y 


x^ 


~  3?  ■ 


siny 


Uc 


sgn 


smy 


m 


~  O'  <  cos  y 


<*OQ  ^ia  sgn  x 


Va^  + 


x^ 


da  >  ~  cos  y  sgn  x 


cos  a  da  0 


sin  a  ,  TT 

da  —  cos  y  sgn  x 


\Ja‘^  — 


x^ 


)  ,  X 


Utilizing  the  previous  results  that  Ji  ~  (1)  const  and  Is  — 0  as  a;  — 0  gives 

^  Q  const  \ 


~  0 


-(1)  ■ 
wy  ^  ——  sm  y  sgn  x 


^  ^  0 


',(1) 


TT 


~  —  COS  y  sgn  x 


(140) 


w 


(1) 


0+ 


=  — TTSini/ 


w 


(1) 


0+ 


=  TT  COS  y 


Thus  the  solution  accounting  for  the  nonlinear  components  generated  at  the  hrst  harmonic 
should  be 

(141) 


^(0)  =  2  ri  (0)  +  Si  (0)  w 


■^a 
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D  Details  of  Linear,  Uniformly  Valid  Solution  to  Two  Layers  of  Uniform  Strat¬ 
ification 

The  governing  equations  for  two  layers  of  uniformly  stratified  flow  in  the  linear  limit  are 


'’P{i),yydx''dx' ,  y  <  H 


'  — oo  J  — oo 


i^{2),yy  + 


'ip{2),yydx'' dx' ,  y  >  H 


'  —oo  J  —oo 


(142) 


'>X(2),x  = 


where  'h  =  ?/  -t-  and  P  =  ^  subject  to 


^(1)  =  -eh{x) ,  1/  =  0 
^(1)  =  ^(2),  y  =  H 
(9^(1)  _  d^(i) 


(143) 


dy  dy 


-,y  =  H 


and  the  radiation  condition. 


Taking  two  derivatives  of  the  equations  for  the  streamfunction  with  respect  to  x  gives 

(^(i),TO  +  +  /^V(i),ra  =  0,  y  <  H 


(^(2), TO  +  ^2^(2))^^  +  /5V(2),to  =  y>  H 


(144) 


and  in  Fourier  space 


-  Pni),yy  =  0i  y  <H 

(^(2), TO  +  ^2^(2))  -  /5^^(2),to  =  y>  H 
^(1),TO  +  Nik^^ii)  =  0,  y  <  H 

(k‘^  -  (3^)  iy2),yy  +  N^k^^2)  =  0,  y>H 

+  (^2!  ^2)^(1)  =  0’  y<H 


(145) 


(146) 


'^{2),yy  + 


Nlk^ 

(A;2  -  /32 


(147) 


-^(2)  =  0,  1/  >  77 
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subject  to  the  boundary  condition  at  the  topography 

^(1)  =  1/  =  0 

and  the  matching  conditions  at  the  tropopause 

^(1)  =  ^(2),  y  =  H 

a^(i)  _  a^(2) 


(148) 


(149) 


dy  dy 

The  general  solution  in  the  Fourier  domain  is 


-,y  =  H 


V’(i)  =  a  (k)  exp  i 


.  Niky 


+  b (k) exp  —i 


.  Niky 


^(2)  =  c  (k)  exp  +  d  (k)  exp  ( -i- 


(150) 


Before  imposing  the  boundary  conditions  at  the  topography  and  tropopause,  it  is  convenient 
to  hrst  take  the  inverse  Fourier  transform  and  impose  the  radiation  condition  in  the  upper 
solution.  In  general,  the  inverse  Fourier  transform  can  be  written 


"/3 


Akx , 


(151) 


because  the  streamfunction  disturbance,  ip,  is  a  physical  quantity  and  must  be  real.  Taking 
the  inverse  Fourier  transform  of  the  general  solution  in  (150)  gives 


^(i)=  23? 


^(2)=  23? 


a  exp 


a  exp  i 


cexp 


'0 


Niky 

y/P^-k\ 

,  Niky  \ 
'  N2ky  ^ 


&exp 


f  -Niky 


AP^-k'\ 


+  6  exp  i 
+  d exp  f 


.  Niky 


y/k‘^  - 
(  -N^ky  V 


N^ky 

I  cexp  I  i- 


'd 


y/k'^  - 


+  d exp  —i 


.  N2ky 


y/k‘^  - 


e^^^dk 


e^^^dk 


e^^^dk 


(152) 


The  radiation  condition  above  the  tropopause  requires  d{k)  =  0  for  k  >  P  and  boundedness 
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as  y  ^  oo  requires  c{k)  =  0  for  k  <  P  giving 


P(i)=  23? 


a  exp 


a  exp  i 


Niky 

^p‘^-k‘\ 

:  N,ky  \ 

Vk'^  -PV 


+  6  exp 


(  -N^ky 


^^p‘^-k‘\ 


+  b exp  —i 


.  Niky 


y/k'^  -  p'^ 


(153) 


,  (  -N^ky 


N2ky 


=  2S  dexp  +  X  "“P  l‘7F^ 


To  facilitate  the  application  of  the  bondary  condition,  we  rewrite  the  solution  in  the  tropo¬ 
sphere  as 


Pii)=  23? 


/  cosh 


/cos 


(  N,ky  ' 
\^P‘^-k\ 

Niky  \  , 


^k^  -  p^ 


+  g  sinh 


g  sin 


Niky 
^p'^  - 

Niky 


e^^^dk 


v/A;2  -  p^ 


e^^^dk 


p{2)  =  23?  <  /  d  exp 


^  ,  (  -N2ky 


e^'^^dk+  I  cexp{i 

h 


WP'^-k\ 


.•  N2ky  \ 


(154) 


y/k^  -  p^ 


e^^^dk 


The  boundary  condition  at  the  surface  of  topography  requires  /  (/c)  =  —eh  {k),  giving 


Pii)=  23? 


-eh  cosh 


-eh  cos 


(  Niky  ' 

[VP^-kp 

Niky  \  , 


g  sinh 


^y/k'^-pp 


g  sin 


Niky 

VP^-k^ 

Niky 


^Pk‘^-pP 


(  -N2ky 


P^^dk 


N2ky 


(155) 


P,2,=  dexp(^^=_le-dA;+  /  cexp  ( 


P'^^dk 
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Then  the  matching  conditions  of  (149)  give 


—eh  cosh 


-eh  cos 


NikH  \ 

NikH 


+  g  sinh 


Ni 


—eh  sinh 


NikH 


g  sin 


+  g  cosh 


Ni  eh  sin 
Solving  gives 
c(h)  = 


NikH  ' 


+  g  cos 


'  NikH  ^ 

(  NikH 

NikH  ' 

V(3^-k\ 

/  NikH  ' 

2Nieh  (k) 


=  dexp  ),  0<k</3 


=  cexp 


y/3‘^-k'\ 

'  iN2kH  ' 


I ,  k  >  (3 


=  -No 


(  -N^kH  ' 


=  iiV2cexp 


'  iNokH  ' 


(N,  -  iV.)  exp  -  (iv.  +  iV,)  exp  (i“S) 

-2Nieh{k) 

(A,.  +  iV,)  exp  (4^^)  +  (iV.  -  iV,)  exp 
W  + /V,)  exp  (MM )  +  (/V.  -  JVO  exp 


d{k)  = 


,  0  <  h  <  /3 

I ,  k  >  j3 

(156) 

{3  <  k  <  oo 

0  <  h  <  /3 


h(/i^)  =  S 


(/V,  +  /V,)  exp  (Ifcgff)  +  (N,  -  JV,)  exp  (=1™^) 

(iV,  -  JV,)  exp  -  (/V,  +  JV,)  exp  (.Mf) 

^  ■  Ve  -  iV.)  exp  (.i^^)  -  (N,  +  N,)  exp  (.^f ) 


eh{k) ,  0  <  k  <  P 


eh{k) ,  P  <  k  <  oo 


Thus  the  streamfunction  solution  is  (155)  with  c(h),  d  {k)  and  g  {k)  given  in  (157). 


(157) 


The  spanwise  velocity  is  determined  from  the  streamfunction  solution.  For  the  linear  case, 
the  spanwise  velocity  is  =  py  Taking  the  Fourier  transform  of  this  gives 


I  f 


(158) 
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Thus  the  spanwise  velocity  can  be  written 


Mj'  (v&)  (t&)l 


—iNi  \  -  /  Niky 


^^2  , _ (  -N2ky 


2*  /„ 


'k'^  -  f3- 


=cexp 


l/>hf 


(159) 
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E  Outer  Solution  for  Two  Layers  of  Uniformly  Stratified  Flow 


The  integrals  for  the  outer  solution  given  in  (98)  are  defined  as 


Jo  =  e3? 


iVcosh  _sinhM^^^ 

Ji  =  e3?  J  /  h  (k) - V 

I  \  ,  /  AThU  \  .  ■,  /  AJUU  \ 


iVcosh  +  sinh 


h{k) 


V  Ufc2-/32  y 


+  i  sin 


iVcosM^  -*sinM^ 


kNjy-H)  \ 


Jg  =  e3?  i  /  h  (fc) 


Jg  =  e3? 


N  exp 


My-H) 

U/32-fc2 


Uo  ^  ^  sinh  +iV  cosh 

I  \y/¥^j 

h  =  e^\  h  (k) - 7 ^  ^  ; - ^e^^^dk 

A  ^  ^  iVcosf  -*sin( 

L  \^k2_^2j  \^k2_p2j 


M  .inh  ( Mkxyo  ^  _  ...h  ( 

-fh  (k)  ^ 

a//52  -  F  iVcoshf  )  +  sinh  f  f  ) 

V  VvA^y  vU/32-fcV 


Jfi  =  e3? 


.AT  sin  f  .  cos  f 

vvfc2-/32y  j  ikx. 


-{h (k)  -^:= - - vv^--y 

^  ^  iVcos[_^j -^sin(_^j 


(160) 


h  =  e3? 


ih  (k^ 


N  exp 


fc(y-g) 

U/32-fc2 


'0  Jp'^  —  k'^  sinh 


TiVcoshH^ 


e^^^dk 


h  = 


h(k) 


AT  exp 

V  Ufc2_^2 


T  Jk^-p^N  cos 


e^^^dk 
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F  Details  -  Unsteady  Wind 

The  kernels  in  the  nonlinear  evolution  equations  are  dehned  as 


1 

Kii  (x,  x')  =  ^  1^  d'^Va  {y'a  +  {y'y'a)i,  -  yy'ai) 


1 

Ki2  (x,  x')  =  ^  dtl^Va  (y'l,  +  (y'y'b)^  -  yyj,^) 


1 

K21  (x,  x')  =  ^  #1/6  (y'a  +  (y'y'a)^  -  yy'ai) 


1 

K22  {x,  x)  =  ^  d^jyb  (y'b  +  {y'y'b)^  - 


(161) 
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